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Abstract 

In the wake of a preceding article 1 30 1 introducing the Schrodinger-Virasoro group, we study its 
affine action on a space of (1 + l)-dimensional Schrodinger operators with time- and space-dependent 
potential V periodic in time. We focus on the subspace corresponding to potentials that are at most 
quadratic in the space coordinate, which is in some sense the natural quantization of the space of 
Hill (Sturm-Liouville) operators on the one-dimensional torus. The orbits in this subspace have finite 
codimension, and their classification by studying the stabilizers can be obtained by extending Kir- 
illov's results on the orbits of the space of Hill operators under the Virasoro group. We then explain 
the connection to the theory of Ermakov-Lewis invariants for time-dependent harmonic oscillators. 
These exact adiabatic invariants behave covariantly under the action of the Schrodinger-Virasoro 
group, which allows a natural classification of the orbits in terms of a monodromy operator on L 2 (R) 
which is closely related to the monodromy matrix for the corresponding Hill operator. 
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Introduction 



The Schrodinger-Virasoro Lie algebra st> was originally introduced in Henkel lfTTl as a natural infinite- 
dimensional extension of the Schrodinger algebra. Recall the latter is defined as the algebra of projective 
Lie symmetries of the free Schrodinger equation in (l+l)-dimensions 

(-2id t -d 2 )iJ/(t,r) = 0. (0.1) 

These act on equation (10.11 ) as the following first-order operators 

Y m = -t m+ W r + i(m+^)t m - l 2r 

M p = \f (0.2) 

with A = 1/4 and n = 0, ±1, m = ±A, p = 0. The Oth-order terms in (10.21 ) correspond at the level 
of the group to the multiplication of the wave function by a phase. To be explicit, the 6-dimensional 
Schrodinger group S acts on ip by the following transformations 

(L_i,L , L x ) : r) -> ifj\t', r) = (ct + d)- [/2 e-? icr2/<ct+d) ip(t, r) (0.3) 

where f = g±£, r' = ^ with ad - be = 1; 

(F ± i) : ifr(t,r) ^ iff(t,r') = e- i{{vt+ro)(l " v/2) il/(t,r) (0.4) 

where r' = r -vt - ro; 

(M ) : <K*,r) ^ e-^ifriU r). (0.5) 

All together these transformations make up a group S, called Schrodinger group, which is isomorphic 
to a semi-direct product of S L(2,R) (corresponding to time-reparametrizations (10.31 )) by a Heisenberg 
group 'Hi (corresponding to the Galilei transformations (10.41) . (10.5I )). Note that the last transformation 
( 10.5b (multiplication by a constant phase) is generated by the commutators of the Galilei transformations 
(10.41) - these do not commute because of the added phase terms, which produce a central extension. 

Now sc - (L n , Y m , M p | n,p e Z, m e A + Z) - made up of all linear combinations of the generators 
corresponding to all possible integer or half-integer indices - is a Lie algebra, as can be checked by 
direct computation. Similarly to the Lie algebra of the Schrodinger group, it is a semi-direct product, 
so - go x b, where go = (L n ) ne z is the centerless Virasoro algebra and I) = (Y m , M p | m € j + Z, p e Z} is 
a two-step nilpotent infinite-dimensional Lie algebra which extends the Heisenberg Lie algebra. It may 
be exponentiated into a group (the Schrodinger-Virasoro group) that we denote by SV. The paper 11301 . by 
C. Roger and the author, studies this Lie algebra for its own sake from different points of view, including 
representation theory, deformations, central extensions. There is a hope that this Lie algebra or related 
ones may help classify strongly anisotropic critical systems and models pertaining to out-of-equilibrium 
statistical physics, notably ageing phenomena, for which the anisotropic dilation (t, r) — > (A 2 t, Ar) (A € R.) 
holds. A systematic investigation of the consequences of a symmetry of the physical system under 
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consideration under the Schrodinger group or related groups has been conducted since the mid-nineties 
(for a short survey, see lfT9l ). 



The starting point for this work is a little different. One of the possible motivations for introducing 
this Lie algebra in the first place is that the group of Lie symmetries of any Schrodinger operator -2id t - 
d 2 + V(t, r) may be represented as a linear combination of the generators introduced in (10.2I ). In other 
words, for any particular Schrodinger operator, the Lie algebra of symmetries is finite-dimensional, but 
the symmetry algebras of all Schrodinger operators are contained in so in the above realization (see 
section 2.5 below for a more precise statement). The proof lies in some sense in a classical paper by U. 
Niederer (see [27]) - who never considered the algebra generated by all possible symmetries. 

Another related way to look at it is that so acts on the space of Schrodinger operators. More precisely, 
SV acts on the affine space of Schrodinger operators with time- and space-dependent potential at most 
quadratic in the space coordinate. We call it S^J '■= {-2id t -d 2 + V2(t)r 2 + V\(i)r + Vo(f)}. It is assumed 
that Vo, V\ and Vj are 2^-periodic in time; this hypothesis is natural when one sets t = e 10 (8 e R/2nZ) 

as a coordinate on the unit circle, so that the generator L n acts as -e m6 de + This restricted space is 

in some sense minimal, which can be seen from the fact that may be expressed in terms of three 
functions of time, just like the elements of so. The phase terms in (10.21) add by commutation with the 
free Schrodinger equation terms of order 1, r and r 2 . One can show that the orbit of any Schrodinger 
operator D e S^f/ has finite-codimension in this space. Hence this space appears to be natural from a 
representation point of view. 

In section 2 below (see section 2.4), we classify the orbits of SV in <S^f/. The classification is 
mainly an extension of Kirillov's results on the classification of the orbits of the space of Hill operators 
under the Virasoro group. These are operators of the type d 2 + u{t). It is well-known (see for instance 
Guieu [13] or [14]) that the group of orientation-preserving diffeomorphisms Diff+(R./27rZ) of the circle 
- which exponentiates the centerless Virasoro algebra - acts on the affine space of Hill operators. Now 
the remarkable fact (despite the apparent differences between the two problems) is that the action of the 
Virasoro group Diff+(R/27rZ) c SV on the quadratic part of the potential, V2(t)r 2 , is equivalent to that 
of Diff+(R/27rZ) on the Hill operator d 2 + V^CO- The reason comes from the fact that the Hill operator 
is the corresponding classical problem in the semi-classical limit (see section 3.2). Hence part of the 
classification may be borrowed directly from the work of Kirillov (see [22]). A. A. Kirillov obtains his 
classification by studying the isotropy algebra Zi<?(Stab„) := {X € L/<?(Diff + (K./27rZ)) | X.(d 2 + u) = 0}. 
There is another equivalent description in terms of the lifted monodromy, which can be explained briefly 
as follows. If (i/q, \p2) is a basis of solutions of the ordinary differential equation (d 2 + u{t))if/{t) = 0, then 
(by Floquet's theory) 



ip 2 {k) + 2n) ) ■ \ tA 2 (f ) 

where M is some matrix (called monodromy matrix) with determinant 1 which does not depend on the 
base point tg. If M is elliptic, i.e. conjugate to a rotation, then the eigenvectors for M are multiplied by 

a phase e l6 . If M is hyperbolic, i.e. conjugate to a Lorentz shift ( 6 _ A j, then the eigenvectors are 

multiplied by a real factor e ±A , hence the solutions of the Hill equation are unstable, going either to zero 
or to infinity when t — > +oo. A nice way to see it (and made rigorous in section 2) is to imagine the vector 
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( \f/ (t) J aS ' rota ^ n §' * n t ' ie Pl ane ma y a l so change norm but never vanishes). The curve described 
by this vector may be lifted to the Riemann surface of the logarithm for instance (obtained from the cut 
plane C \ R_), so that it turns by an angle unambiguously defined in R. This gives the lifted monodromy. 

The space S a J/ has been considered independently by mathematicians and physicists, with similar 
motivations but different methods (that turn out to be equivalent in the end). The general idea was to 
solve the evolution problem associated with D e S^/, i.e.to show that the Cauchy problem Dtfr = 
with initial condition tA(0, r) = (AoW has a unique solution and compute it explicitly. The usual method 
in mathematical physics for such time-dependent problems is to consider the adiabatic approximation: 
if one puts formally a small coefficient e in front of d t , the problem is equivalent by dilating the time 
coordinate to the equation (-2i<9 f - d 2 r + V(st, r))ij/ = 0, so that V is a potential that is slowly varying 
in time. Suppose that A e (?) := -d 2 + V(st, r) has a pure point spectrum [A n (t),n e N} for every t, 
where A n is C°° in t, say, and let i// n (t) be a normalized eigenfunction of A £ (t) satisfying the gauge- 
fixing condition (i{/ n (t),ip n (t)) = 0. Then there exists a parallel transport operator W(s,t) carrying the 
eigenspace with eigenvalue A„(s) to the eigenspace with eigenvalue A n (t), and a phase operator ®(s, t), 
given simply by the multiplication by a phase e3 / Ani - s) ds on each eigenspace, such that the solution of 
the Schrodinger equation is given at first order in s by the composition of W and <!>. One may see the 
solutions formally as flat sections for a connection (called Berry connection) related in simple terms to the 
phase operator (see [4j). This scheme may be iterated, giving approximate solutions to the Schrodinger 
equation that are correct to any order in e (see for instance Joye l20l ). but it is rarely the case that one 
can give exact solutions. By considering the related classical problem, G. Hagedorn (see ifToTl ) constructs 
a set of raising and lowering operators (generalizing those associated to the usual harmonic oscillator) 
for general Schrodinger operators in S^/, and uses them to solve the equation explicitly. The same 
set of operators had been considered previously by two quantum physicists, H. R. Lewis and W. B. 
Riesenfeld (see |[25l ). and obtained by looking for an exact invariant, i.e. for a time-dependent operator 
/(?) (not including the time-derivative) such that % = % + i[I(t), \(d 2 - V(t, r))] = 0. They find for each 
operator D in S?JJ a family of invariants (called sometimes the Ermakov-Lewis invariants, see ||28Tl ) 
depending on an arbitrary real solution £ of a certain differential equation of order 3 (see Proposition 
3.1.4), constructed out of generalized raising and lowering operators and spectrally equivalent to the 
standard harmonic oscillator -k(d 2 - r 2 ). These invariants have been used to solve quite a few physical 
problems, ranging from quantum mechanics for charged particles to cosmology (see ifTTll . lfl2l . ll28l . 
|[29l for instance). It turns out that very few Schrodinger operators have an exact invariant of the type 
I(t) = f2(t, r)d 2 + f\(t, r)d r + fo(t, r). These may be expressed, as shown by H. Lewis and P. Leach (see 
O), in terms of three arbitrary functions of time (the exact expression is complicated). Exact invariants 
allow in principle to solve explicitly the original problem, at least if one knows how to diagonalize them 
(which is the case here). Hence (provided one requires that an exact invariant exists) the space S"/^ is 
maximal. 

There are three new features here: 

- the action of the Schrddinger-Virasoro group on S a J^ (which is essentially a conjugate action, leav- 
ing all invariant quantities unchanged, for instance the spectrum and the monodromy) makes it possible 
to reduce the study to five families of operators, with qualitatively different properties (see section 2.4). 
They are mainly characterized by the monodromy of the associated Hill operator d 2 + V<2(?), but there 
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also appeal - some non-generic orbits in cases when the quadratic and linear parts of the potential are 
'resonant'. The non-periodic case is much simpler, since (locally in time) all Schrodinger operators in 
«S<2 w& formally equivalent (see section 3 below). The coefficients of the Ermakov-Lewis invariants are 
related in a very simple way to the invariants of the orbits; 

- one is interested in Schrodinger operators with time-periodic potential. Hence one may consider (as 
in the case of ordinary differential operators, see above) the monodromy, which is a bounded operator 
acting on L 2 (R). The monodromy operator is given explicitly and shown to be closely related to the 
classical monodromy of the related Hill operator; 

- the computation of the monodromy in the case when the associated Hill operator is hyperbolic(see. 
above) requires the use of an Ermakov- Lewis invariant associated to a purely imaginary function 
which is equivalent to the standard harmonic 'repulsor' -^(d 2 +r 2 ). The reason (explained more precisely 
in section 3 below) is that the usual Ermakov-Lewis invariants are defined only if /y 2 (£) > 0, where the 
invariant quantity /y 2 (£) (quadratic in £) is associated to the Hill operator d 2 + Vi(t) and its stabilizer %(t) 
in Lie(Diff + (R/27rZ)). The stabilizer satisfies a linear differential equation of order 3 and has generically 
only one periodic solution (up to a constant). If one does not require £ to be periodic, then /„(£) may be 
chosen to be positive, which is perfectly suitable for a local study (in time) but is of little practical use 
for the computation of the monodromy. If however one requires that £ be periodic, then is negative 
in the hyperbolic case, unless one chooses £ to be purely imaginary. Hence one is naturally led to use 
the spectral decomposition of the harmonic 'repulsor' (which has an absolutely continuous spectrum 
equal to the whole real line). Usually there is no adiabatic scheme, hence no phase operator, in the case 
when eigenvalues are not separated by a gap. But in this very particular case, such a phase operator 
may be computed and is very analogous to that obtained in the elliptic case, for which the spectrum is 
discrete. There exists also some non-generic cases (corresponding to a unipotent monodromy matrix for 
the underlying Hill operator) for which 7y 2 (£) = 0. The natural invariant is then spectrally equivalent 
either to the bare Laplacian -\d 2 or to the Airy operator -\(d 2 - r). 

One of the main results may be stated as follows (see sections 3.4,3.5,3.6): the monodromy operator 
is unitarily equivalent to the unitary multiplication operator f(k) — » e lkT ~ my f(k), where y is some con- 
stant and k is the spectral parameter of the model operator -\(d 2 r + Kr 2 ) (k — ±1, 0) or -\(d 2 - r), and 

T - 4jk reai ) or i Jq K 155 imaginary). The above integrals must be understood in a general- 
ized sense if £ has some zeros; a complex deformation of contour is needed then. Comparing with the 

usual Berry phase e~- X Ak< -^ ds , one sees that the eigenvalue X^t) - -2k is constant, but that the natural 
(possibly singular) time-scale is r := J 

The paper is organized as follows. 

Section 1 is preliminary and contains notations and results (contained in QUI ) concerning the Schrodinger- 
Virasoro group and its action on Schrodinger operators. 

Section 2 is dedicated to the classification of the orbits and of the isotropy subgroups Go '■= {g £ 
SV | g.D = D), D e S^fJ (see section 2.4). It contains long but necessary preliminaries on the action 
of the Virasoro group on Hill operators. The connection to the results of U. Niederer is made in the last 
paragraph. 

We solve the monodromy problem for the Schrodinger operators of the form -2i<9 ; -d 2 + V2(t)r 2 + y 
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(y constant) in section 3. We study first the corresponding classical problem given by the associated Hill 
operator, x + V%(t)x = (an ordinary differential equation). The solution of the quantum problem is then 
easily deduced from that of the classical problem. In either case, the monodromy is obtained by relating 
the Ermakov-Lewis invariants to the orbit data. 

Finally, we show in section 4 how to parametrize a general Schrodinger operator -2id t -d^ + Vi{t)r 1 + 
V\(t)r + Vo(0 e *S<2 by means of a three-dimensional invariant (g(t), 5\(t), (52(0) (see Definition 4.2). 
The parametrization is one-to-one or 'almost' one-to-one depending on the orbit class of the potential V2 
(S V-orbits in S?fJ have generically codimension 2, whereas adjoint orbits corresponding to the invariant 

have generically codimension 2 or 3). The action of the Schrodinger- Virasoro group on S a J 2 ^, once 
written in terms of the invariant, becomes much simpler, and is easily shown to be Hamiltonian for a 
natural symplectic structure. A generalized Ermakov-Lewis invariant may also be written in terms of 
this three-dimensional invariant. We then solve the monodromy cases for the 'resonant' cases left from 
section 3. 

Notation: The notation Vect(S ! ) := {^(t)d t | £ € C°°(]R/27rZ)} will be used for the Lie algebra of C°°- 
vector fields on the torus R/2nZ. The infinite-dimensional group Diff + (R/27rZ) of orientation-preserving 
diffeomorphisms of the torus R/2nZ (also called centerless Virasoro group) has a Lie structure, and its 
Lie algebra is Yect(S l ) (see lfT4l for details). 

Let us gather here (for the convenience of the reader) a few notations scattered in the text. Time and 
space coordinates are usually (at least starting from section 2) denoted by 9 and x (see explanations before 
Lemma 1.6 for the passage to Laurent coordinates (t,r)). Stabilizers in Vect^ 1 ) of the Hill operator 
dg+ Vi(0) are usually denoted by £ (which is either real or purely imaginary). If £ is purely imaginary, then 
one sets £ := in. As for (operator) invariants of the Schrodinger operators (see section 4), we write them 
as \ [a(#)x 2 - b(9)d 2 x - ic(6)(xd x + d x x) + d(6)(-id x ) + e(ff)x + . The correspondence between the 
vector invariant (£,61,62) and the operator (generalized Ermakov-Lewis) invariant is given in Theorem 
4.4. 

1 The Schrodinger- Virasoro group and its action on Schrodinger opera- 
tors 

We recall in this preliminary section the properties of the Schrodinger group proved in ||3"U1 that will be 
needed throughout the article. 

Definition 1.1 (see ED, Definition 1.2) 

We denote by 5v(k), k - or \, the Lie algebra with generators L n , Y m , M n (n e Z, m e k + Z) and 
following relations (where n, p e Z,m,m' e k + Z) : 

[L n ,L p ] = (n-p)L n+p 
n 

[L n , Y m ] = (- - m)Y n+m , [L n , Mp\ = -pM n+p ; 
[Y m , Y m >] = (m- m')M m+m >, 
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[Y m ,Mp]=0, [M n ,Mp\ = 0. 

Iff (resp. g, h) is a Laurent series, f = £„ eZ /„/' I+ \ msp. g = Z nei(+ z8nt n+ K h = EneZ V> then 
we shall write 

£ f = J^f n L n , 2/ g = J^g n Y n , M h = Yj h n M "- ( u ) 

Note the shift in the indices in the Laurent series which disappears in the Fourier coordinates, see 
remarks preceding Lemma 1.6 below. 

It is often unimportant (or a matter of taste) in this paper whether the shift index k is or I (see 
remarks after Theorem 2.4.2 though). In this section (unless otherwise stated) so stands indifferently for 
st>(0) or so(l/2). In the following sections, we shall abbreviate so(0) to so for convenience. 

Definition 1.2 (see (301, Definition 1.3) 

Denote by dn A the representation of so as differential operators of order one on K. 2 with coordinates 
t, r defined by 

i i 

dn A (L n ) - -f +l d t - -(n + \)frd r + -i(n + l)«?" _1 r 2 - {n + \)Af 

dn A (Y m ) = -t m+ W r + i(m + r 

dn x {M p ) = it? (1.2) 

Proposition 1.3 (see Il30l , Theorem 1.1) 

1. The Lie algebra st> can be exponentiated to a Lie group denoted by S V. It is isomorphic to a semi- 
direct product SV - Go x H, where Gq - Diff + (R/27rZ) is the group of orientation-preserving 
diffeomorphisms of the torus R/27rZ, and H - C°°(R/27rZ) X C°°(]R/27rZ) (as a vector space) is the 
product of two copies of the space of infinitely differentiable functions on the circle, with its group 
structure modified as follows: 

(a 2 ,p 2 )-(a l ,pi) - (ax + a 2 ,P\ + f3 2 + -(a[a 2 - a ia ' 2 )). (1.3) 

The semi-direct product is given by: 

(1; O) = (0; (1-4) 

and 

(<f>; 0).(1; (a,P)) = (0; ((<*') s(a o <p),p o c/>)). (1.5) 

2. The infinitesimal representation dn A ofsv can be exponentiated to the following representation of 
the group SV on C°° functions of two variables, 
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(i) 

(^(0;O)/)(fV) = (0'(*))-V7co r /(*,r) 
jf € Diff + (R/27rZ) induces the coordinate change (t, r) — » (?', /) = (0(0, r $'(!)); 

(ii) 

{n A (\;{a,P))f){t'y) = ^i^^,) 

if(a,fS) e C°°(]R/27rZ) X C°°(]R/27rZ) induces the coordinate change (t,r) — > (?, r') - (t,r - 
a(t)). 

Definition 1.4 (see ED, Definition 2.1) 

Le? «S' m fte vector space of second order operators on R 2 defined by 

D e <S'" 7 <=> D = /j(-2i^ - d 2 ) + V(f, r), fc, V e C°°(R 2 ) 

and S a ^ c «S' m f/ze o^wie subspace of 'Schrddinger operators' given by the hyperplane h — 1. 

/ti of/zer wonis, arc element of S a ^ is the sum of the free Schrddinger operator —2id t — <9 2 aw<i o/a 
potential V. 

Proposition 1.5 (see 11301 . Proposition 2.5, Proposition 2.6) 

Le? cr 1/4 : 5 V — > Hom(S hn ,S bn ) the representation of the group of SV on the space of Schrddinger 
operators defined by the left-and-right action 

<ri/4fe) : D -» ^5/4 {g)Dn i / 4 (g) - ' , g e S V, D € 

r/ie« o"i/4 restricts to an affine action on the affine subspace S a ^ which is given by the following formu- 
las: 

<n/ 4 (0;O).(-2icf ; -d 2 + V(t,r)) = 

-2id t -dj + <p'(t)Vm),r^[f(fj) + K 2 ®(cf>)(t) (1.6) 

0-1/4(1; {a, b)).(-2id t -d 2 r + V{t, r)) = 
-2\d t -d 2 .+ V(t, r - a(t)) - 2ra"{t) - {2b' (t) - a{t)a"{t)). (1.7) 

where ^7 — 5 ("I 7 ) * J tne Schwarzian derivative. 

One may also consider a generalized left-and-right action cr A (g) : D — > 7r / i + i(^)D7r /i (g)" 1 , but then 
the subspace «Sj (see Definition 2.1.2) is not preserved by cr /i | D iff + ( R / 2 ^z) any more, which ruins all 
subsequent computations. Actually 1 /4 corresponds to the 'scaling dimension' of the Schrodingerian 
field in one dimension (see ED). 



We shall occasionally use the time-reparametrization 

<p : R/2;rZ -» S 1 a tf(l), -> f - e ie (1.8) 
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from the torus to the unit circle. It allows to switch from the Fourier coordinate 9 to the Laurent coordi- 
nate t. In particular, 

£ t „+t = tti/4(0;O) £ e inB 7ri/4(0;O) _1 , J/ 1 = 7Ti/ 4 (0; 0) J/ e i«e tti/ 4 (0; 0) _1 , = tti/ 4 (0; 0)ALme tti/ 4 (^ 

(1.9) 

If « is an integer, J/ + i should be understood to be acting on the two-fold covering of the complex 
plane where the square-root is defined; conversely, if n is a half-integer, then }f e u,9 acts on 47r-periodic 
functions. In other words, the 'natural' choice for so should be 50(5), resp. sr>(0) in the Laurent, resp. 
Fourier coordinates. 

Applying formally the formulas of Proposition 1.3, one gets 

(ffiMfcOr 1 /)^,*) = (k" i9 ) 1/4 e"^/(e ie , ±xe i( ^) (1.10) 

(with some ambiguity in the sign) which is an 87r-periodic function. Applying now (still formally) 
Proposition 1.5 yields the following result, which can be checked by direct computation. 

Lemma 1.6 

Let f(t, r) be a solution of the Schrddinger equation 

{-2id t -d 2 r + V(t,r))f(t,r) = 0. 

Then 

f : (0, x) -» e~ m e-* ixl f(e i(> , xe m2+!T/4) ) (1.11) 
is a solution of the transformed Schrddinger equation 

-2id e - d 2 x + \x 2 + ie ie V(e ie , xe im+7T/4) ) f(9, x) = 0. (1.12) 



In the following sections, we shall (except when explicitly mentioned) always work with the Lie 
algebra st>(0) in the Fourier coordinates 9, x (i.e. the Lie algebra generated by the X/, }f s and Wt/, with 
In-periodic functions f,g,h), and write so instead of so(0) for simplicity. 

2 Classification of the Schrddinger operators in S // 

From now on, whe shall concentrate on the affine subspace of Schrodinger operator with potentials which 
are at most quadratic in the space coordinate. As mentioned in the Introduction, this subspace is invariant 
under the action of SV. The purpose of this section is to classify the orbits. 

2.1 Statement of the problem and connection with the classification of Hill operators 

Let us first define two natural subspaces of S a ^ '. 

Definition 2.1.1 (Schrodinger operators with at most quadratic potential) (see El, Prop. 2.6) 
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LetS a J 2 f = {-2id g -dl+V2(6)x 2 + V l (6)x+V (6)} c S aff be the affine space of Schrddinger operators 
with a potential which is In-periodic in time and at most quadratic in the coordinate x. 

Definition 2.1.2 (Schrddinger operators with quadratic potential) 

Let S a ^ = {-2idg - d 2 + V2{9)x 2 } c S a J/ be the affine space of Schrddinger operators in S"^ with 
time-periodic potential proportional to x 2 . 

We do not assume V2 to be positive. Hence what we really consider are harmonic 'oscillators- 
repulsors', corresponding to the quantization of a classical oscillator-repulsor with time-dependent Hamil- 
tonian \(p 2 + V2(9)x 2 + V\(6)x + Vq{9)). If V\ = 0, then t he classical equation of motion = 
-V2(6)x - \V\{9) has as an attractive, resp. repulsive fixed point depending on the sign of V2. If V2 
is not of constant sign, things can be complicated; it is not clear a priori whether solutions are stable or 
unstable. We shall come back to this problem (which turns out to be more or less equivalent to the a 
priori harder quantum problem, at least as far as monodromy in concerned) in section 3.2. 

The first subspace is preserved by the action of S V (see Proposition 1 .5) and is in some sense 
minimal (the SV-orbit of the free Schrodinger equation, or of the standard harmonic oscillator -2i3<9 - 
d 2 + a 2 x 2 , contains 'almost' all potentials which are at most quadratic in x). As we shall prove below, 
the orbits in S a ^ have finite codimension. 

Let us write down for the convenience of the reader the restriction of the action of 0-1/4 to S a J^: let 
D = -2id e -d 2 x + V 2 (9)x 2 + Vi(0)x + V (6), then 

<ri /4 (# 0)(D) - -2id e -d 2 x + {^' 2 . V 2 o <p + ^0(0)1 x 2 + (<f>' 3/2 . V l o <f>) x + </>' . V o <f> (2. 1) 

fff f/ ^ 

- recall 0(0) = ^- - \ (^-J is the Schwarzian derivative -, and 

o- i/4 (l; {a, b)){D) = -2id e -d 2 x + V 2 x 2 + (Vi - 2aV 2 - 2a")x + {Vq - aV\ + a 2 V 2 - 2b' + aa"), (2.2) 
while the infinitesimal action is given by 

do- 1/4 (£ f )(D) = -(If" + 2fV 2 + fV' 2 )x 2 - (fV[ + \f'Vi)x - (fV' Q + fV ), (2.3) 
do- m & g + M h )(D) = -2(g" + gV 2 )x - (2h' + gVi). (2.4) 

These four formulas are fundamental for most computations below, and we shall constantly refer to 
them. 

Similarly, S a ^ is preserved by the cri/4-action of Diff + (R/27rZ) (see Proposition 1.5). It turns out 
that the orbit theory for this space is equivalent to that of the Hill operators under the Virasoro group. Let 
us first give some notations and recall basic facts concerning Hill operators. 

Definition 2.1.3 

A Hill operator is a Sturm-Liouville operator on the one-dimensional torus, i.e. a second-order 
operator of the form oi + u(G) where u(9) e C°°(]R/27rZ) is a 2n-periodic function. 
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The action of the group of time-reparametrizations on a Hill operator may be constructed as follows. 
Starting 'naively' from the simple action of diffeomorphisms on functions, 

ip^ipocp, (f>€ Diff+(R/27rZ), 

one sees that (<9 2 + u)(ip) = is equivalent to the transformed equation (d 2 + p(9)d + q(9))(ip o cp) = if 
one sets p = and q = cp' 2 . u o (p. Then one uses the following: 

Definition 2.1.4 (Wilczinsky's semi-canonical form) (see Magnus-Winkler, l26l . 3.1, or Guieu, 
US, Proposition 2.1.1) 

Ifipisa solution of the second-order equation (d 2 + p(6)d + q(9))ip = 0, then tp := A(9)ip is a solution 
of the Hill equation (d 2 + u(9))ip = provided 

A(9) = exp|^ ^ p{s) dsj (2.5) 

for some % and 

1 , 1 2 

u = --p - -p +q. (2.6) 

One obtains in this case A = (cf>')~ 1 ^ 2 , and the transformed operator reads: d 2 + (cf>') 2 . u o cp + i©(0), 
where is the Schwarzian derivative. The presence of this last term shows that this transformation 
defines a projective action of Diff+(R/27rZ). Summarizing, one obtains: 

Proposition 2.1.5 (see Guieu, [ 13 ] or Guieu-Roger, [ 14]) 

The transformation 

d 2 + u — > (p t (d 2 + u) := d 2 + (cp') 2 . u o <p + i©(0) 

defines an action q/"Diff + (R/2ffZ) on the space of Hill operators, which is equivalent to the affine coad- 
joint action on vix\ (i.e. with central charge c = ^). A solution of the transformed equation may be 

2 

obtained pom a solution tp of the initial equation (d + u)ip = by setting (p t ip = (<p')~zip o <p. In other 
words, the solutions of the Hill equations behave as (-^-densities. 

The important remark now is the following: 
Lemma 2.1.6 

The above action o/Diff + (R/27rZ) on the space of Hill operators is equivalent to the cr\j^-action of 
Diff + (R/2;rZ) on the space S a 2 ff . 

Namely, Proposition 1.5 above (see also (12.11 )) shows that 

o- l/4 (4>){-2id -d 2 x + V 2 (6)x 2 ) = -2id e - d\ + V 2 (9)x 2 , (2.7) 

where the potential V2 is the image of V2 (viewed as the potential of a Hill operator in the coordinate 9) 
by the diffeomorphism <p, i.e. (f>*(dg + V%(ff)) = d 2 e + VziO)- Once again, this should not come as a surprise 
since the Hill equation is the semi-classical limit of the Schrodinger operator (see section 3.2). □ 
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So we shall need to recall briefly the classification of the orbits of Hill operators under the Virasoro 
group. There are mainly three a priori different classifications, which of course turn out in the end 
to be equivalent: the first one is by the lifted monodromy of the solutions (see for instance B. Khesin 
and R. Wendt, [21]); the second one consists in looking for normal forms for the solutions, either an 
exponential form for non- vanishing solutions or a standard form for a dynamical system associated with 
the repartition of the zeros (see the article by V. F. Lazutkin and T. F. Pankratova, E31 ): the third one, 
due to A. A. Kirillov (see 11221 ) proceeds in a more indirect way by looking at the isotropy groups. We 
shall need the first and the last classification for our purposes. They are the subject of the two upcoming 
subsections (see also for a related review and application to the global Liouville equation).. 

2.2 Classification of Hill operators by the lifted monodromy 

Let us now turn to the classification of the orbits under the Virasoro group of the space of Hill operators. 

Consider a pair (t/q , ^2) of linearly independent solutions of the Hill equation (d 2 + u)\p = 0. It is 
a classical result (a particular case of Floquet's theory for Schrodinger equations with (space)-periodic 
potential) that 



for a certain matrix M{u) € 5L(2, R) (independent of 9), called the monodromy matrix. Starting from 



a different basis ~ , one obtains a conjugate matrix M(u). The above action of the Virasoro group 



on the Hill equation leaves the monodromy matrix unchanged, as can be seen from the transformed 
solutions 0*^1,0*^2- Hence the conjugacy class of the monodromy matrix is an invariant of the Hill 
operator under the action of the diffeomorphism group. 

Floquet's theory, together with the orbit theory for S L(2, R), imply that d 2 + u is stable (meaning that 
all solutions are bounded) if |TrM| < 2 or equivalently, if M is elliptic, i.e. conjugate to a rotation matrix; 
unstable (meaning that all solutions are unbounded) if |TrM| > 2 or equivalently, if M is hyperbolic, i.e. 

/ e A \ 

conjugate to a Lorentz shift I _^ I, A > 0. If \TrM\ = 2, then M can be shown to be conjugate 

/ 1 27T \ , 

either to ±Id or to the unipotent matrix ± I ^ I; in the latter case, + u is semi-stable, with stable 

and unstable solutions. Two linearly independent 27r- or 47r-periodic solutions exist when M = ±Id; only 
one in the unipotent case; and none in in the remaining cases. 

An important result due to Lazutkin-Pankratova (see [23]) states that all stable Hill operators are 
conjugate by a suitable time-reparametrization to a Hill operator with constant potential d 2 + a, a > 
0. They also distinguish between oscillating and non-oscillating equations (oscillating equations have 
solutions with infinitely many zeros, while non-oscillating equations have solutions with at most one 
zero), but we shall not need to go further into this. Let us just remark that (as they also show) non- 
oscillating operators are also conjugate to a Hill operator with constant potential d 2 + a, with a < this 
time. Hence operators of type II, resp. Ill of Kirillov's classification (see Definition 2.3.4 below) are 
exactly the unstable, resp. semi-stable oscillating operators. 




(2.8) 
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A complete classification of the orbits under the action of Diff+(R/27r/Z) may be obtained by consid- 
ering the lifted monodromy. Set I J^j^ J = Af(w)(0) I ^® J . The path 9 -> M(w)(0) e 5L(2, R) may 

be lifted uniquely to a path 9 -> M(u)(9) e SL(2, R) such that M(w)(0) = Id, where SL(2, R) is the univer- 
sal covering of 5L(2, R.). This procedure defines a unique lifted monodromy matrix M(u) := M(u)(2n) 
modulo conjugacy. 

The following arguments (see ETI ") show briefly why this invariant suffices to characterize the orbit 

of u under diffeomorphisms. Set I "f 1 ^ |= J£(6) [ C ° S0J ^ |. The Wronskian 

\ $z(Q) I \ sinw(9) / 

W := ^1^2 - 

(a constant of motion) is equal to <x/(0)£(0), hence to' - j- is of constant sign, say > (by choosing 
W > 0). By the action of Diff(R/2;rZ), one can arrange that a>' is constant, while a>(0) and (o(2n) 

remain related by the homographic action of M(u), viz. cotana>(27r) = ccotan^( ( 0)+rf if M ^ ~ ( ^ ^ j e 
5L(2,R). The lifting of the monodromy produces a supplementary invariant: the winding number n := 
L(co(2/r) - w(0))/2^] = ^ ^1 (L ■ l=entire part), namely, the integer number of complete rotations 
made by the angle oj. 

This change of function is particularly relevant in the elliptic case. Choose a basis I 1 ) such that 



M = ( C ° S ] ~ Sin / V Then ±^ = W (2;r) - cu(0) = W jjf C&r] - 
\ sini cosi / J() w 

/ e A \ 

If M = I _ j I is hyperbolic instead, set rather 



<A2 



tftfi) = \w)\e 2a>{6 \ = ^(0)l<T 2 " (9) (2.9) 

with f(0) = 2(tffnff 2 )(ff), so that ±i = w(2ff) - w(0) [2«r]. Then one finds oj' = — |, hence w = -W J 

The functions | and a; are not well-defined if \p\ or i/r2 has some zeros. Supposing u is analytic, the 
functions ^1,^2 may be extended analytically to some strip O = {|Im0| < e}. Choose some contour 
T c Q. avoiding the zeros of \fj\ and ip2 such that (assuming £(0) + 0, otherwise use a translation) 
T(0) = and T(27r) = 27r. The idea is to keep T real away from some symmetric neighbourhood U £ of 
the zeros, and to complete the path with half-circles centered on the real axis of radius s around each 
zero, taken indifferently in the upper- or lower-half plane (compare with section 3.2 below where more 
care is needed). Suppose >I/\{9q) = for instance, so (/^(^o) - a ^ and ^2(^0) = Then 

C d9 C e ° +£ d9 0~0 O C e ° +E d6 n 
-W =-W - = -Wp.v.\ ±1- (2.10) 

Jme - E M 0+£] m X- e 0-00 + io m J d0 - £ m 2 

(depending on the position of the half-circle with respect to the real axis) since g_^ +i0 = P-v.^zj" + 
(see for instance flU) and the residues of ^ at the zeros of £ are It is clear from the above 
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definitions that £ has only simple zeros, in even number. Hence —W J r ^ = -W p.v. ^ = A [in]. 
By exponentiating, one obtains a monodromy matrix in PS L(2, R) - SL(2, R)/{±1}. 

Finally, if M is unipotent, M = ± [ n ^ ) in some basis ( ^ ), set ifri(9) = 0)1^2(6) and £ - so 



1/ \ ^2 

that a»(27r) - a»(0) + a. Then a/ = so oj is once again defined as - W J -Mr if £ does not have any 
zero. In the contrary case, one uses a deformation of contour as in the hyperbolic case, to obtain 

C d6 f e ° +£ d6 (<9-<9 ) 2 
-W \ T7^ = -W — — — 7 - : ■ (2.11) 



rn[e - £ ,eo +£ ] ^) X-s (0-9o±W ty\(ff) 

Since (e^TIoF = P- V -WW ± in5 8o and T^f = 1 + " 9q)2) ~ Since ^' ^ = -ViiBoWiiOo) = -, 
the Dirac term does not make any contribution at all this time, hence 

a = co(2n) - co(0) = -W J^||, (2.12) 

where Y : [0, 2n\ — > C is an arbitrary contour as defined above. 

Summarizing: 

Proposition 2.2.1 (see ED for (ii)) 
(i) The lifted monodromy of the operator d 2 +u is characterized by the (correctly normalized) quantity 

Cm^IrWr where f eStab - 

f ii) The orbits under the diffeomorphism group of the space of Hill operators are characterized by 
the conjugacy class of their lifted monodromy. More precisely, the lifted monodromy defines a 
bijection from the set of orbits onto the space of conjugacy classes of{SL(2,R)\{±l})/{±l}(an 

element M e S L(2, R) has to be identified with its opposite -M). 
2.3 Kirillov's classification of Hill operators by isotropy subgroups 

Another classification, also useful for our purposes (and more explicit in some sense), is due to Kirillov. 
Introduce first 

Definition 2.3.1 

LetStab u , u € C 00 '(R./2ttZ) be the isotropy subgroup (or stabilizer) ofd 2 +u in Diff + (R/2ffZ), namely, 
Stab u := [<p € Diff + (R/27rZ) | $ t (d 2 + u) = d 2 + u). (2.13) 



Proposition 2.3.2 (definition of the first integral /) (see |fT3l ) 
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1. 



Let £; € C°°(K./27rZ): then ^ € Lie(Stab u ) if and only if ^ satisfies 

l -r + 2< + viz = o. 



(2.14) 



2. Let I u (^) := — \<;' 2 + 2w£ 2 . Then I u (^) is a constant of motion if £ e Lie(S tab u ). 

3. Consider <p € Diff+(R./27rZ) and the transformed potential u such that (f>*(d 2 + u) — d 2 + u. Then 

h(<p'- 1 -W) = (2-15) 

4. Consider the Hill equation (d 2 + u)ip(8) - 0. If(ft\ , 1//2) is a basis of solutions of this equation, then 
£ := a\irf/ 2 + IxxYityxtyi + anft\ (a\\,a\2> a 22 e R) satisfies the equation 

" + 2ug + u't = (2.16) 

In other terms, £ e Stab u is in the isotropy subgroup of the Hill operator d 2 + u. 
Conversely, any solution of (I2.16i can be obtained in this way. 

5. (same notations) consider in particular £, = ip 2 + ip 2 ^. Then /„(£) = W ifW is the Wronskian of 
(ft i, ft2), namely, W - ft\ft' 2 - ft[ft2 (constant of the motion). 

Note (see 3.) that (ft' 1 ■ £ ftjd is the conjugate of the vector £d e Vect(S l ) by the diffeomorphism 
ft. Hence one may say that the first integral / is invariant under the (adjoint-and-coadjoint) action of 
Diff + (R/2/rZ). 

Consider now the (adjoint) orbit of £ under Diff+(R/2nZ). Clearly, Jjj ^ (if well-defined, i.e. if £ 

has no zero) does not depend on the choice of the point on the orbit since T^r^f^^ - ^ ■ It 
is easy to see from Prop. 2.3.2(2) that £ either never vanishes (case I), or has an even number of simple 
zeros (case II), or has a finite number of double zeros (case III). Cases II, III correspond to a hyperbolic, 
resp. unipotent monodromy matrix (see discussion in section 2.2). In case II, /„(£) = - j£'(to) 2 < if 

to is any zero. The principal value integral p.v. j Q $z is well-defined. In case III, /„(£) = and the 
regularized integral J r -M? (see above) is well-defined and independent of the choice of the contour T. 
Note that A. Kirillov uses instead the following regularization, ]im s ->o f . v ^ — § (where U £ is a 
symmetric e-neighbourhood of the zeros) with C chosen so that the limit is finite. The two regularizations 
are different. Both are perfectly satisfactory to define an invariant of the orbits, but computations show 
that the Berry phase is proportional to J j^. 

Now the integral -M: (case I) and its variants for case II, III are invariants under the diffeomor- 
phism group. The discussion in section 2.2 shows that they characterize the lifted monodromy of d 2 + u. 
The invariant is also needed to fix u uniquely in case I (see Prop. 2.3.2(2)) since £ stabilizes all 
operators of the type d 2 + u + p- (C e R). It turns out that - or its variants - and I u (0 (in cases II and 
III), together with a discrete invariant neN, suffice to distinguish between the different adjoint orbits 
of stabilizers (note that general adjoint orbits may be much more complicated, see Ifl4ll ). One has the 
following: 

Proposition 2.3.3 - Classification of the coadjoint invariants and of the orbits (see Kirillov |[22l ) 
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1. Case I: £ is conjugate by a diffeomorphism <f> to a ( non-zero ) constant adg, a + 0. Hence <f>' 1 . £ o 
(p € Lie(S tab i)i +a ) for a certain constant a. The stabilizer Stab d i +a is: 

— («) 

(i) (non-generic case) either isomorphic to SL (2,R) (the n-fold covering of S L(2,W)), with 

2 

Lie(S tab$2 +a ) = RSgOR cos n8dg(BM. sin n8dg if a = \for some n e N*; then the monodromy 
in PSL(2,R) = SL(2, R)/{ + l} is trivial, while the lifted monodromy matrix is the central 
element in SL(2, W)/{±1] corresponding to a rotation of an angle nn; 

(ii) or (generic case) one-dimensional, equal to the rotation group Rot c Diff + (R/27rZ) gener- 
ated by the constant field d in the remaining cases. 

The invariants are given by - 2aa 2 , - The monodromy can be in any 

conjugacy class ofPSL(2, R) except ±Id. 

2. Case II: g is conjugate to the field a sin n9( 1 + a sin n9)dg, n = 1 , 2 . . ., < a < 1, which stabilizes 
d + u n a , where 



n 

Un,a(9) '■= 



1 r 1 + 6a sin nG + Aa 2 sin 2 n6 



(1 + asin?i0) 2 



(2.17) 



The monodromy matrix is hyperbolic. The invariants take the values I u (^) - —2a n < 0, 

3. Case III: g is conjugate to ^± iH , a '■= ±(1 + sinn^)(l + asinn6)d, < a < 1, corresponding to a 
potential v n ^ a , 

" (or - 1 ) 2 + 2a(3 - a) sin n9 + 4a 2 sin 2 n6 



v n ,a(0) 



n 2 



(2.18) 



(1 + a sin n9) 2 

The monodromy matrix is unipotent. The invariant I u (g) vanishes, while f j^L: = ~ 2n . The 

J i s+W (l-a)Vl-ff 2 

discrete invariant n suffices to characterize the orbit of d 2 + u. 

In cases II and III (provided a > 0), the stabilizer is one-dimensional, generated by £dg. 

In the generic cases (case I, a + n 2 /4, n = 0, 1, ... or case II) the monodromy matrix is elliptic, 
resp. hyperbolic, if and only ifl u (g) > 0, resp. /„(£) < 0. In cases I (a = 0) and III (with unipotent 
monodromy), = 0. 



There is a mistake in Lemma 3 of 11221 (the potential w n>a given there is not correct). The potential 
v n>a was missing, together with the value of J r -fj^- Both are obtained by straightforward computations. 

This classification is also natural when one thinks of the behaviour of the solutions (see Lazutkin- 
Pankratova ||23l and section 2.2). In particular, case II (resp. Ill) correspond to operators with unstable 
(resp. semi-stable), oscillating solutions, while case I corresponds to operators with stable, oscillating 
solutions (a > 0), resp. unstable, non-oscillating solutions (a < 0), resp. semi-stable, non-oscillating 
solutions (a = 0). 

Note that in the case I generic, the three-dimensional isotropy subalgebra contains fields £ of type I, 
II (a = 0) and III (a = 0), hence the following nomenclature: 

Definition 2.3.4 
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Ifd+u has a stabilizer i; of type I, or of type II, III with a = 0, then d 2 + u may be turned into a Hill 
operator with constant potential, and we shall say that the operator d 2 + u (or the potential u) is of type 

I. If d 2 + u has a stabilizer of type II, resp. Ill with a ± 0, then we shall say that d 2 + u and u are of type 

II, resp. type III. 

Similarly, we shall say that the Schrodinger operator —2idg — d 2 + V2(6)x 2 + V\(9)x + Vq(6) is of type 
I (resp. II, III) if the Hill operator d 2 d + V2(6) is of the corresponding type. 

Note that the cases I generic (a + '-j-, n = 0, 1, . . .) and II are generic (i.e. dense in S^). 



Now the eigenvalues of the monodromy matrix (and also the lifted monodromy) can easily be ob- 
tained once one knows the values of the invariants ^ and /«(£). The following Lemma gives the 
link between the two classifications: 

Lemma 2.3.5 

Suppose D = d 2 + u is of type I ( with a ± 0) or II (i.e. its monodromy is either elliptic or hyperbolic ). 
If D is of type I non generic, conjugate to d 2 + n 2 14 for some n > 1, choose £, to be conjugate to some 
non-zero multiple ofdg. Now (in all cases) normalize g by requiring that I u (^) = 1, so that £; is real in 
the elliptic case and purely imaginary in the hyperbolic case. Then the eigenvalues of the monodromy 
matrix are given by exp ±i J Q or eX p + 1 p V j Q 

Proof. 

Coming back to the discussion in section 2.2, one checks easily (with the normalization chosen 
there) that /„(£) = 2W 2 in the elliptic case, and /«(£) = -2W 2 in the hyperbolic case. Choose a basis of 
solutions (^1,^2) sucn that W = 1 and multiply £ by i in the hyperbolic case. Then (in both cases) the 
eigenvalues of the monodromy matrix (±iA in the elliptic case, and ±A in the hyperbolic case) are given 
by exp ±i ^ or the exponential of the corresponding principal value integral. □ 



2.4 Classification of the SV-orbits in S~? 2 J 

This problem can be solved by extending the above results, which may be interpreted as the decompo- 
sition of S"f^ into Diff + (]R/27rZ)-orbits. Let us first compute the stabilizers of some operators that will 
be shown later to be representatives of all the orbits. We choose to present the results in the Fourier 
coordinates (9, x). The orbits of type I, resp. Ill split into orbits of type (i), (i)bis, resp. (iii), (iii)bis due 
to the presence of the linear term V\(9)x in the potential. 

The computations depend on the formulas of Proposition 1.5, see formulas d2.il) , ( 12.21 ), ( 12.31 ), ( 12.41) 
for more convenience. 

Definition 2.4.1 

If D € «S< { , we denote by Go the stabilizer of D in the Schrddinger-Virasoro group SV, i.e. Go - 
{geSY\ o- m (g).D = D). 

Recall the notation Stab„, u e C°°(]R/27rZ) is used for the stabilizer in Diff + (R/27rZ) of the corre- 
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sponding Hill operator. 

Note that Mi = Mo (whose exponential amounts to the multiplication of the wave functions t// by a 
constant phase) acts trivially on any operator D, hence M\ e Go always. The rotation group 8 — > 9 + Gq 
generated by do-\j$(£]) = dcri/4(Lo) - -de will be denoted by Rot. 

In the following classification, we shall call harmonic oscillators (resp. harmonic repulsors) opera- 
tors with elliptic, resp. hyperbolic monodromy. 

(i) Time-independent harmonic oscillators or repulsors 

Set D a ^ y := -2ide-dl+ax 2 +y (a, y e R). It is clear that L_i -de leaves D a ^ 7 invariant in all cases. 
Suppose first for simplicity that y - 0. Then Go - (Go)z)»*#z) (see Proposition 1.3 for notations) is 
a semi-direct product, so one retrieves Kirillov's results (see Proposition 2.3.3, case I) for (Gq)d', 
to be specific, Lie((G )D n2/40 ) = R<5<? © R-£ s i n „e © RXcosne if n e N\ and Lie((G () ) Da0 ) = Rd e 
otherwise. 

Now (1; (a, b)) e Hjj a0 if and only if b' - and a" - -aa. The latter equation has a non-trivial 
solution if and only if a = (in which case Lie(Ho) - RJ/i © RAti) or a = n 2 /4, n > 1 with « 
even, in which case Lie(Ho) = R2/ C osne/2 © RJ^sin«6»/2 © -Mi. Then exp ^Xi c 5 , L ( " ) (2, R) is the 
rotation of angle 2n, while exp ^adXil[LiefK ),LieCK D )] is a rotation of angle n. 

The isotropy groups Go are the same in the case y t 0, except for a different embedding in- 
volving sometimes complicated components in the nilpotent part of SV which do not change the 
commutation relations (so that Go is no more a semi-direct product (Gq)d « Hp). 

All together, one has proved: 

Theorem 2.4.2 

1. If a = n 2 /4, where n > 2 is an even integer, then Gd ii2/4 - S l} \2, R) x 'Hi is isomorphic to 

— (n) 

an n-covering of the Schrddinger group; the semi-direct action of SL (2, R) quotients out 
into an action of the two-fold covering S L (2, R). The Lie algebra of the group S L (2,R) 
acts as Rde © R(X s inn6» + M_i ysinng ) © R(X cos «6» + ^-i yC osn6»)- A f ter transformation to the 
Laurent coordinates (t, r) (and supposing y = 0), Go 2/40 is the connected Lie group with Lie 
algebra (L Q ,L ±n ) x (Y ±n/2 ,M Q ) c st>(0). 

2. If a = n 2 /4, where n > 1 is odd, then Gd m2/40 - Sl (n) (2, R) x exp RMi. 

3. If a = 0, then G Df)y = Rot x exp(R}/i © RAti) - (R/2/rZ) x R x (R/2ttZ) is the commutative 
group of constant translations-phases. After transformation to the Laurent coordinates (t, r), 
it is the connected Lie group with Lie algebra (Lq, Yq, Mq) c sd(0). 

4. In the generic case a + n 2 /4, n = 0, 1, . . . one has simply Gd = Rot X exp RMi - (R/2ttZ) 2 . 

It is natural in view of these results to consider the two-fold covering H {2) of H obtained by 
considering 4/r-periodic fields. Then the stabilizer in S V := Go x # (2) of D n 2/ 4 (n > 1 odd) is 

isomorphic to S L (n) (2, R)k-Hi as in the case of an even index n. This time Lie(!H\) - (Y ±n /2, Mo) c 
SD (i). 
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The best-known case is a - 1/4 (n - 1), y - 0. In the Laurent coordinates (t, r), writes 
-2id t - d 2 , namely, it is the free Schrodinger equation. Then S L(2, R) x 'Hi acts on D1/4 in the 
usual way (see formulas d0.2| |) in the Laurent coordinates. 



(i)bis Special time-independent harmonic oscillators with added resonant oscillating drift 



Consider 



D = -2ide - d\ + n 2 x 2 + C cos(n8 - <r/2) . x + y 



(C,o-,y e R, C ± 0, n > 1 integer). Then computations show that G D =* RxRx R/2nZ is 
three-dimensional, generated by 

-Cl-cos(2«6l-0-) + Mc. sin3(n6»-0-/2) ~~ M , c OS (2n0-cr) t y cos ( 2 „0_ (r ))' (2-19) 

^ sM +^^ s(M _ ff) (2-20) 

and M\. One checks (by direct computation) that the value of the associated invariant I n 2(l - 
cos(2n# - cr)) is 0. 

(ii) Time-dependent Ince harmonic repulsors of type II 
Consider 

D„, a ,y = -2id e - d 2 x + u n ^ a (8)x 2 + y, n = 1,2, . . . , a e (0, 1) (2.21) 



where 

1 + 6a sin n8 + 4a 2 sin 2 n6 

u n,a(9) — T 



(1 + a sin nd) 2 



(2.22) 



Then (see preceding subsection) - € Lie(Go) (£ ^ 0) if and only if £ is proportional to 
£ n ^ a , with f„ >Q: - sin«0(l + asinn8)dg. Now 



^0-1/4(0//, + M f2 ).D = 

if and only if / 7 ' - and /" + u iua f\ = 0. The latter equation is known under the name of Ince 's 
equation (see Magnus and Winkler, [26 ]). The change of variable and function 8 — > 5(8) = f -n| , 
-> y(^) = (1 + o-cos 25) fc/4ff /i (£(#)) with b = -2a[l + -4=] turns the above equation into 
the standard form 

( 1 + a cos 25)y" + b sin 28y' + (c + d cos 26)y = 

2 2 9- 

with a = a, c = 1 - 7^-7, d = a\3 + t^-? + -p=l. Conditions for the coexistence of two 
independent periodic solutions of Ince's equation have been studied in detail. In our case, there is 
no periodic solution since d 2 + u, ha is unstable (see discussion in section 2.2). Hence 

G D = exp(R(X^ - |A%,.„) 6 RMi) - R x (R./2/rZ). (2.23) 

(hi) Non-resonant time-dependent Schrodinger operators of type III 
Consider 

D n ,a, y = 2id e -d 2 x + v niCt (8)x 2 + y 
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(n - 1,2,,.., a € (0,1)). See equation Q2.18I ). Similarly to case (ii), - e Lie(Go) 

(£ + 0) if and only if £ is proportional to %±, n ,a, where £±, n ,a = ±(1 + sin«0)(l + a sinnG)dg. Then 
rfo"i/4(J/y, + Mf 2 ).D = if and only if fL = and f" + v n>a f\ - 0. This is once again Ince's 
equation, with parameters a - a, b - -2a, c — 1 - 2a, d - 3a. One verifies immediately that 
y{8) = cos 5 is the unique (up to a constant) periodic solution of this semi-stable Hill equation, 
corresponding to f\{6) := (1 + asmnG)^ cos(| - Note that f+ >njQI = /j 2 (so that /i is - up to a 
sign - the unique C°° square-root of £+, n ,a)- Hence 

G D = exp (rU^ ~ \M{ ±Jia )©Ry /l ©RM 1 )-]RxRx (R/2;rZ) (2.24) 

(iii)bis Schrodinger operators of type III with added resonant drift 
Consider 

D = -2idg - d 2 x + v n ^(G)x 2 + C(l + asinn9) l/2 cos(- - n-)x + y 

(C ^ 0) with v^ a as in case (iii). Set £(8) = (1+sin nG)(\+a sin nQ) and f{6) = {l+a sin nG)^ cos(7r/4- 
n6/2). Recall £ = / 2 . 

Suppose Xf + + At/ 2 stabilizes D. Then (see ( 12.31 )) 

2(/T + v^/D = C(£/' + ^'/) - 4C/ 2 /' (2-25) 

and 

/a = + chf)< ( 2 - 26 ) 

The kernel of the operator d 2 + v n>a is one-dimensional, generated by /. Hence the above equation 
(|2T25T > has a solution if and only if jf^/' + \g f)f d9 = 0, which is true since (£f + - 
(/ )'. Now equation (12.261 ) has a solution if and only if f\ is chosen to be the unique solution 
orthogonal to the kernel of d 2 + v„ jQ ., namely, if j Q f\f dQ = 0. 

Now 

d(Ti/ 4 (}/ gl + M g2 ).D = -2(g" + v„, a gi)x - fgi - 2g' 2 (2.27) 

vanishes if and only if g\ - f (up to a multiplicative constant) and J Q gif dG - 0. The two 
conditions are clearly incompatible. 

All together one has proved that Go = exp (R(£ ( i +S j n „ex l +a sin nS) + M j\ + Mf 2 ) © RMij - Ix 
R/2jtZ (with /i,/2 solving equations (12.25 1 ). (12.261 )) is commutative two-dimensional. 

Explicit but cumbersome formulas for fufz w& easv to derive from the proof of Lemma 4.10 
below. We shall not need them. 

There remains to prove that we have classified all the orbits in S^fJ. 
Theorem 2.4.3 

Any Schrodinger operator D in S"/^ belongs to the orbit of one of the above operators. 
Proof. 
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Let D € <S a J 2 . Suppose first that V2 is of type I. Then one may assume (by a time-reparametrization) 
that V2 = a is a constant. The operator D belongs to the orbit of D a>y (case (i)) for some y if and only 

.(2) 

if V\ - 2{a" + aa). If a + n (or rr/4 if one considers the S V -orbits) then this equation has a unique 
solution for every Vi. If a = n 2 , then a Fourier series V\ = Yjk cue lke is in the image of d 2 e + a if and only 
if c±n — 0- This analysis accounts for the two cases (i), (i)bis. 

Suppose now V2 is of type II. By a time-reparametrization one may choose V 2 = u na . The operator 
D belongs to the orbit of D, ua ^ (see case (ii)) for some 7, provided V\ - 2{a" + u n>a a). Since d 2 + u n ^ a 
(acting on C°°(]R/27rZ)) has a trivial kernel, it has a bounded inverse and the unique solution of the above 
equation is C°°. Hence D belongs to the orbit of D n ^ y . 

Finally, suppose V2 is of type III. One is led to solve the equation V\ = 2{a" + v n ^ a a). Recall 
- (1 + a sinn0) 1,/2 cos(? - n|) solves the equation /" + v n ^ a f\ = 0. Hence V\ - 2(a" + v n ^ a a) 

has a solution if and only if L n Vi(9)(\ + a sin n6) l/2 cos(| - n|) d6 - 0, which accounts for cases (iii), 
(iii)bis. □ 

Note that Schrodinger operators of type III are generically of type (iii)bis, and Schrodinger operators 
of type I with a = n 2 , n = 1, 2, . . . are generically of type (i)bis. 

Corollary 2.4.4 

2 

For generic orbits ( type ( i) with a + n > 0, or type ( ii) ), the isotropy group is two-dimensional, 



given by exp R(£ f + 3J Sl + Mg 2 ) © RMi - Ix R/2ttZ or R/2nZ x R/2nZfor some triple (^, 81 , 8 2 ) with 
£ t 0. 



Let us finish with a remark. Consider a potential V2{0)x 2 + V\(6)x + Vo(^) of type (i), (ii) or (iii). As 
we shall see in the next section, the monodromy of the corresponding Schrodinger operator depends only 
on the (conjugacy class of the) invariant £ and the value of the constant y (which acts as a simple energy 
shift). Computing the invariant £ is a difficult task in general, but suppose it can be achieved. How does 
one determine the constant y ? We give an answer for generic elliptic or hyperbolic potentials of type (i). 

Lemma 2.4.5 

Let D - —2ide — d 2 + V2(9)x 2 + V\(6)x + Vo(6) be of type (i), elliptic or hyperbolic, generic, so that 
D is conjugate to a unique operator D a ^ y - —2ide — d 2 + ax 2 + y (a € R, a + n = 0, 1, . . .). Then y 
may be retrieved from 



Start from the model operator D a ^, with stabilizer £ = 1, and apply successively 0-1/4(0; (0, 0)) and 
o"i/4(l; (g, h)). Then one obtains the operator 




(2.28) 



where W\ is the unique solution of the equation (d 2 + V2)W\ = V\. 



Proof. 



D := -2id e - d\ + V 2 (9)x 2 - 2 ((d 2 + V 2 )g) x + (y4> + g{d 2 + V 2 )g - 2h) 
(see formulas (IZTTl.dZIl)). Now * 0(0) dO - 2n since e Diff+(R/2nZ). Hence the result. 



□ 
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2.5 Connection to U. Niederer's results 



We are refering to a classical paper by U. Niederer (see |27]) concerning the maximal groups of Lie 
symmetries of Schrodinger equations with arbitrary potentials. One may rephrase his main result as 
follows (though the Schrodinger- Virasoro had not been introduced at that time). U. Niederer shows that 
any transformation 

if/(f,r) -> M*,r) = expif g (g~\t,r))ip(g~\t,r)) 

(where g : (t, r) — > (t' , r') is an arbitrary coordinate transformation and f g an arbitrary 'companion func- 
tion' corresponding to a projective action) carrying the space of solutions of the Schrodinger equation 

(~2id t - d 2 + V(t, r))ifr(t, r) = (2.29) 

into itself in necessarily of the form a\ /4(g) for some g e SV. This is - by the way - an elegant way 
of introducing the SV group in the first place. Then Niederer gives a necessary and sufficient condition 
for g € cri/4(SV) to leave (12.291) invariant, and produces some physically interesting examples. Let 
us analyze some of these examples from our point of view. It should be understood that Niederer's 
examples are given in the Laurent coordinates (?, r) and should hence be transformed by using Lemma 
1.6 to compare with our results. 

(i) V - (free Schrodinger equation): this case corresponds after the transformation at the end of 
Section 1 to the potential V(9, x) = |x 2 , with invariance under the full Schrodinger group (see case 
(i) in subsection 2.4, with a = 1/4 and y = 0). 

(ii) V = -gr (free fall) corresponds to V(9, x) - \x 2 - ge l(e ^ 2+3n / 4 ^ x (a 47r-periodic potential), which 
belongs to the same orbit as case (i) (free Schrodinger equation in the Laurent coordinates). 

(iii) V = \uj L r 1 (harmonic oscillator) may be obtained from the free Schrodinger equation by the time 
reparametrization t(u) = tan cou for which the Schwarzian derivative is a constant, &(t) = 2oj 2 (see 
formulas in Proposition 1.5). 

(iv) V = k/r 2 (inverse-square potential), corresponding to the operator -2i<9 e - d 2 + ^ + kx~ 2 (har- 
monic oscillator with added inverse-square potential) in the Fourier coordinate. The operator is 
not in S^f, but the (time-independent) inverse-square potential is interesting in that this is the 
only potential left invariant by all transformations V(t, r) — > (p'{t)V{(p{t), r ^<p'{t)) (see formulas in 
Proposition 1.5). So this equation is invariant by the kernel of the Schwarzian derivative, i.e. by 
the homographic transformations. 

3 Monodromy of time-dependent Schrodinger operators of non-resonant 
types and Ermakov-Lewis invariants 

We 'solve' in this section all Schrodinger operators in S^J of class (i), (ii) or (iii) by using the Ermakov- 
Lewis invariants, to be introduced below. Since any such operator is conjugate to an operator of the type 
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-2idg - d 2 + V2(0)x 2 + y (y constant), and y corresponds to a simple energy shift, we shall implicitly 
assume that the potential is simply quadratic (V\ = Vo = 0). 

Lemma 3.1.2 and Proposition 3.1.4 yield explicitly an evolution operator U(6,Gq), i.e. a unitary 
operator on L 2 (R) which gives the evolution of the solutions of the Schrodinger equation from time 6q 
to time 6. This operator gives the unique solution to the Cauchy problem and allows to compute the 
(exact) Berry phase. The arguments in Lemmas 3.1.1, 3.1.2 and Proposition 3.1.4 are reproduced from 
the article of Lewis and Riesenfeld ([25 ]). Unfortunately this method gives the monodromy only in the 
elliptic case (i.e. for operators of class (i) with a > 0). So we generalize their invariants to the hyperbolic 
and unipotent case; the invariant we must choose in order to be able to compute the monodromy is not 
a harmonic oscillator any more, but an operator with absolutely continuous spectrum. Nevertheless, 
it turns out that there does exist a phase operator, given in terms of the (possibly regularized) integral 
X ~W) ^ or a certam stabilizer £ of the quadratic part of the potential. The key point in order to get the 
whole picture is to make the bridge between Kirillov's results and the Ermakov-Lewis invariants. 

3.1 Ermakov-Lewis invariants and Schrodinger- Virasoro invariance 

Let H - i(— d 2 + V2(9)x 2 ) be the (quantum) Hamiltonian corresponding to a time-dependent harmonic 
oscillator. The evolution of the wave function if/(G, x) is given by: ideip(Q, x) - H\p(G, x), or Dif/ = 
where D = -2id e + 2H = -2id e -d 2 x + V 2 (G)x 2 . 

The Ermakov-Lewis dynamical invariants were invented in order to find the solutions of the above 
equation. The idea is simple. Suppose 1(6, x) is a time-dependent hermitian operator of the form 
Y^ = Qlj{G,x)d x which is an invariant of the motion, i.e. j f I = d t I + i[I,H] = 0. Suppose also that, 
for every fixed value of 6, 1(6, x) (defined on an appropriately defined dense subspace of L 2 (R, dx), for 
instance on the space of test functions) is essentially self-adjoint and has a purely point spectrum. For 
simplicity, we shall assume that all multiplicities are one, and that one may choose normalized eigenvec- 
tors which depend regularly on 0, namely, 

7(6>, x)h n (9, x) = A n (9)h„(d, x) (3.1) 

and f \h n (Q, x)\ 2 dx = 1. The fact that I is an invariant of the motion implies by definition that Iif/ is a 
solution of the Schrodinger equation if iff is. The following lemma shows how to solve the Schrodinger 
equation by means of the invariant I: 

Lemma 3.1.1 (see l|25l ) 

1. The eigenvalues A n (6) are constants, i.e. they do not depend on time. 

2. \in + m, then (h m (6), (id e - H)h n (9)) = 0. 

Proof. 

(i) Applying the invariance property |^ + \[I(9), H(6)] -Oto the eigenvector h n (G) yields 

dl 1 

—h n (G) + -(1(G) - A n (6))H(9)h n (6) = 0. 
oG l 
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Taking the scalar product with h m (9) gives a first equation, 

(h m {9), ^-h n (6)) + r(/U0) - A n (8))(h m (9), H(9)h n (9)) = 0. (3.2) 
39 1 

The eigenvalue equation I(9)h n (9) = A n (9)h n (9) gives after time differentiation a second equation, 
namely 

—h n (6) + (1(6) - A n (9))h n (9) = A n (9)h n (9). (3.3) 
06 

Combining equations ( 13.21 ) and (13.31 ) for n - m yields A n (9) - 0. 
(ii) Combining this time equations (13.21 ) and (13.31 ) for n + m yields the desired equality. 

The above Lemma shows that one may choose eigenvectors h n (9) that satisfy the Schrodinger equa- 
tion by multiplying them by an appropriate time-dependent phase, which is the content of the following 
Lemma. 

Lemma 3.1.2 

Let, for each n, a n (9) be a solution of the equation 

^ = (h n (9), (id e - H)h n (9)). (3.4) 
d9 

Then the gauge-transformed eigenvectors for the invariant I 

~h n (9) = e ia " (6) h n (9) (3.5) 

are solutions of the Schrodinger equation. 

In other words, the general solution of the Schrodinger equation is: 

m := J] c n e iaM h n (9) (3.6) 

n 

where c„ are constant (time-independent) coefficients. 

Let us specialize to the case when H is a time-dependent harmonic oscillator as above, i.e. H = 



\(-d 2 x + V2(9)x 2 ). A natural idea is to assume the following Ansatz 



1(9) = i [-b(9)d 2 x + a(9)x 2 - ic(9)(xd x + d x x)] . 



2 

This problem has a unique family of non-trivial solutions: 
Definition 3.1.3 (Pinney-Milne equation) 

The non-linear equation 



t + f(9X-^=0 (3.7) 
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(K > 0) is called a Pinney '-Milne equation. If K - 1, then we shall say that (13.7\) is a normalized 
Pinney-Milne equation. 

Of course, every Pinney-Milne equation can easily be normalized by multiplying the function by the 
constant factor K l/4 . 



The following Proposition summarizes results due to Lewis and Riesenfeld (see [25"]). 
Proposition 3.1.4 (Ermakov-Lewis invariants for time-dependent harmonic oscillators) 



1. The second-order operator &£,(£ 2 ) 



S£(t)(8) = - 



(3.8) 



is an invariant of the time-dependent harmonic oscillator —2ide — d 2 x + V2(8)x 2 provided £ is a 
solution of the following normalized Pinney-Milne equation : 



£+v 2 (0)<r-l = o. 



(3.9) 



Setting g = £ 2 , one may also write equivalently 



i 

2? 



X 2 + (tfd x + -i;x) 2 

£ 2 i 

-fd 2 x + (1 + ^-)x 2 + -tf(xd x + d x x) 



(3.10) 



2. Set 



and 



a{9) := — 

V2U(0) 



a* (9) = — 

(formal adjoint of the operator a(9)). Then 



- (ae)d x + W) X ) 



+ (md x +iao)x) 



&C(£>(0) = a*(9)a(9) + -. 



(3-11) 



In other words, for every fixed value of 9, the operators a(9), a* (9) play the roles of an annihilation, 
resp. creation operator for the (time-dependent) harmonic oscillator &£.(£). 

3. The normalized ground state of the operator a(9) is 



1 1 



1 1 



(3.12) 
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4. The solutions of equation (13.4\) giving the phase evolution of the solutions of the Schrddinger 
equation are given by 

1 dff 

a„(6) = -(n + -) J — (3.13) 
provided one chooses the time-evolution of the eigenstates h n by setting 

(h n ,d e h n )= 1 -{n + ^){a-?). (3.14) 

The above choice for the time-evolution of the eigenstates appears natural if one requires the standard 
lowering and raising relations a(0)h n (6) = mh n -\{9), a*{6)h n {9) - (n + l)2/i n+1 (#). Then computations 
show that 

(h n ,d e h n ) = (h ,d e h ) + i^{U-l 2 ). (3.15) 

Hence there only remains to choose the time-evolution of the ground-state ho. This particular choice 
leads to the (n + ^-factor typical of the spectrum of the harmonic oscillator. Note that the h„(6) do not 
satisfy the gauge-fixing condition typical of the adiabatic approximation (see A. Joye [20] for instance). 
But this phase choice leads to a nice interpretation of the phases a n (up to a constant) as a canonical coor- 
dinate conjugate to the classical invariant £X C / (see Lemma 3.2.2 below) for the corresponding classical 
problem, in the generalized symplectic formalism for which time is a coordinate, so that the problem be- 
comes autonomous (see Lewis-Riesenfeld ESI ; see also section 4 for the symplectic formalism). Also, 
as mentioned in the introduction, the natural time-scale (both for the classical and the quantum problem) 

The connection with the preceding sections is given by the following classical lemma (see [26 ], chap. 
3), which is an easy corollary of Proposition 2.3.2: 

Lemma 3.1.5 



1. Let ^ be a (non-necessarily periodic) solution of the equation 

if" + 2< + ^ = 0, (3.16) 
so that ^ stabilizes d 2 + u. Then £ := V? is a solution of the Pinney-Milne equation 

£' + u{GX - y = (3. 17) 

where /«(£) ; - — j£' 2 + 2u£ 2 is the constant defined in Prop. 2.3.2 (2). 

In particular, if = ifn+if/L where (iff i, iff 2) is a basis of solutions of the Hill equation (d 2 + u)ij/ = 0, 



and £ - V?> then 

W 2 
1 

where W := iff 1 iff 2 — iff [iff 2 is the Wronskian of the two solutions. 



C' + u(9X- — =0 (3.18) 
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2. Consider £ € Stab u such that £ - satisfies the Pinney -Milne equation \3.18\ , and a time- 
reparametrization (p. Then £- := <p'~ l . £ o (p is a stabilizer of d 2 + u := (p*(d 2 + u) and X := 
satisfies the transformed Pinney-Milne equation 



~ W 2 
C + K - = 



/or same constant W. 



be the associated Ermakov-Lewis invariant. 



The interesting point now is that one can choose the Ermakov-Lewis invariant in such a way that the 
invariant associated to the image of the time-dependent harmonic oscillator D by a time reparametriza- 
tion (through the representation 0-1/4) is its image by a very natural transformation (essentially, by the 
corresponding change of coordinates). This provides an elegant, natural explanation for the complicated- 
looking phase appearing in the formulas for 0*1/4. 

Theorem 3.1.6 

Let D := —lidg -d 2 + V2(0)x 2 be a time-dependent harmonic oscillator, g satisfy the Pinney equation 
C' + V 2 {-j, = 0, and&m 2 ) = \ [(f) 2 + {%d x + Cx) 2 

Let <p € Diff + (R/27rZ) be a time-reparametrization and V2 be the image of V 2 through <p, defined by 
0-1/4(01? - -2id e -d 2 + V 2 {9)x 2 . 

Then: 

1. X := ((p' o <p~ l )~- . £ o <pT x satisfies the transformed Pinney equation 1" + V2I — 4$ = 0. 

2. Consider the transformed Ermakov-Lewis invariant 

1 



££(f )(x) := - 



l) +{ild- x + 'ix) 2 

a vt de 



(3.19) 



where (9, x) - ((p(9), x ^(p'{9)) are the transformed coordinates. 
Then 

££(f ) = 7ri / 4(0££(^ 2 )7r 1/4 (0)- 1 . (3.20) 
In particular, ££(£ 2 ) is an Ermakov-Lewis invariant for o~\/4((p)D. 

Proof. 

1. follows from Lemma 3.1.5 (2). This implies that ££.(l 2 ) is an Ermakov-Lewis invariant for 
o"i/4(0) • D. Supposing one has proved that S£(( 2 ) is the conjugate of ££(£ 2 ) by ^1/4(0), then it follows 
once again that B£,(X 2 ) is an invariant for 0-1/4(0) . D since 

(<ri/4(0 . D)S£(? )-&l{l 2 ) (0-1/4(0) .D) = <f>' 7ri /4 (0)D .S£(?) .n m {<py l -n m {cp) .S£(?) . cp'Dn^Y 1 = 
(the function of time cp' commutes with the operator &£(£ 2 )). 
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So all there remains to show is that ££(<T 2 ) is indeed conjugate to &£(£ 2 ). This is actually true for 
both terms appearing inside parentheses in the expression for the Ermakov-Lewis invariant (and trivial 
for the first one). Set £ = i£d x + £'x and £ = \t,d x + ^ Jc. Then a simple computation shows that 

, 1 d>" 
£ = \Cd x + xt + -x—r 

On the other hand, 

1 ma U <l>"(0) , .2 I ... \: <P"V) ,2 , \ 

- -ms)T l/A e^-m x me)d x + {'(fix) [fitfUe-VTm* ^{9\ x ^W))) • (3.21) 
Hence 7Ti/4(0)£7ri/4(0)~* — £. □ 



We now want to be able to write the general solution of the Schrodinger equation as 

>A(0) = J e iak(ff) c k h k (6)d(r(k) 

(for some spectral measure cr on a set S, a discrete measure in the case studied by Lewis and Riesenfeld) 
with periodic eigenstates h k and a phase a k with periodic derivative, i.e. given by integrating a periodic 
function, so that 

iff(6 + 2n) = J e iAk e iak{e) c k h k {6)do-{k) 

where the A k :- a k (9+2n)-a k (6) are constants and measure the rotation of the eigenstates h k after a time 
In. Then the monodromy operator is unitarily equivalent to the multiplication operator f(k) — > f(k)e lAk 
on L 2 (Z, do~). 

Consider any Schrodinger operator with quadratic potential V2(6)x 2 and an associated non-zero vec- 
tor field £ e Stabiy-i) as before. (We postpone the discussion of 'resonant' operators (classes (i)bis 
and (iii)bis) to the next section.) It turns out that the eigenstates h k and the measure a can be taken 
as the (possibly generalized) eigenfunctions and spectral measure of one of the three following 'model' 
operators H, depending on the sign of the invariant /«(£): 

(i) (4(£) > 0) : take for H the standard harmonic oscillator 

H = -^(d 2 x -a 2 x 2 ) (aeW); 

this case corresponds to harmonic oscillators of type (i), i.e. Schrodinger operators of type (i) 
conjugate to -2idg - d 2 x + a 2 x 2 with a 2 > 0; 
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(ii) (I u (£) < 0) : take for H the 'standard harmonic repulsor' 

H = -j(d 2 x + a 2 x 2 ) (aeR); 

this case corresponds to harmonic repulsors of type (i), i.e. operators of type (i) conjugate to 
-2ide - d 2 x - a 2 x 2 (-a 2 < 0), and operators of type (ii); 

(iii) (/ M (£) = 0) : take for H the usual one-dimensional Laplacian, 

this case corresponds to operators of type (i) conjugate to the free Schrodinger operator -2idg-d 2 , 
and operators of type (iii). 

Note that this classification is equivalent to the classification of the (conjugacy classes of) mon- 
odromy matrices for the associated Hill operators d 2 e + V2(9) into elliptic, hyperbolic and unipotent ele- 
ments. 

The next section circumvents the spectral analysis technicalities by solving the associated classical 
problem. The essentials for understanding the (operator-valued) monodromy for the quantum problem 
are already contained in the study of the (S L(2, R)-valued) monodromy of the ordinary differential equa- 
tion x = -V2(6)x, so we found this short digression convenient for the reader. Then we study the spectral 
decomposition of the above model operators. Finally, we solve the quantum problem for a quadratic 
potential V2(0)x 2 and compute the monodromy operator. The general case D e S"^ may be reduced to 
the quadratic case D e after applying some transformation in SV, except for the operators of type 
(i)bis and (iii)bis which will be treated in the last section. 



3.2 Solution of the associated classical problem 

The associated classical problem (obtained for instance as the lowest-order term in % in the usual semi- 
classical expansion) is a Hill equation. 

Definition 3.2.1 (classical problem) 

Let H be the classical hamiltonian H = \{p 2 + V2(6)x 2 ). 

The asociated motion in phase space reads x - d p H - p, p - -d x H - -V2, which is equivalent to 
the Hill equation (dj + V 2 )x(8) = 0. 

Lemma 3.2.2 



1. Suppose V2 is of type I with a + or of type II, and choose g e StabV^ so that /«(£) = 2 is real 
in the elliptic case and purely imaginary in the hyperbolic case). Then 



&ZcM){x) ■- l - 



^ + ^ x - l J_ x) 2 



(3.22) 



is an invariant of the motion. 
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2. Suppose V2 is of type I with a = or of type III ( so that the associated monodromy is unipotent), 
and take any ^ e StabV^ % £ 0. Then 



ax- --x) 2 

2£ 



(3.23) 



is an invariant of the motion. 



Proof. Simple computation (&£ c i may be obtained from the quantum Ermakov-Lewis invariant by 
letting h go to zero). □ 

Assuming V% is elliptic, i.e. of type I with a > 0, one may choose £ > 0. Then the equation 

£-£c/(£)W = C, C constant is equivalent to ( j&j + z 2 = C after the function- and time-change t(6) = 

j 6 ||py, x{9) = ^(9)z(r(9)), with obvious solutions cos r, sin t. Hence a basis of solutions of the equation 
of motion is given by 

1 r® do' 1 r® d9' 

*i(0) = ^(0) cos J _ * 2 (0)=£5(0)sin J — (3.24) 

r ■ . i^N . . CO riff Cd Hff rr,. / X l \ ,r. n / COS ^ ~~ Sm T 

Assume for instance that £(0) = 0, and choose J - J Then I \{2n) - 



X2 J \ sin T cos T J 

with 7 = Jq* Hence the eigenvalues of the monodromy matrix are given by ±iT . 

In the hyperbolic case (type I with a < 0, or type II), £ := irj is purely imaginary. The above formulas 
(13.241 ) give solutions of the Hill equation on either side of any zero of £ (note that the normalization 
/„(£) = 2 implies £(#) ~e->e ±2i(9 - 6q) near any zero, so that ( 13.241 ) defines a continuous function, as 
should be of course), but the easiest way to define the solutions X\,X2 globally is to use a deformation of 
contour. One may always assume that £ is analytic on some complex neighbourhood of R (it is conjugate 
by a time-reparametrization to some u n>a which is entire, see Prop. 2.2.3). Define a contour T from to 
2jt which avoids the zeros of £ by going around them along half-circles of small radii centered on the real 
axis. This time (see discussion in section 2.2), the half-circles must be chosen alternatively in the upper- 

and lower-half planes so that Re £(z) > on T. Then ( %x ) (2tt) = I ° OS "L Sm „T ) • I Xl 1(0) as 



X2 J \ sin T cos T ) \ X2 

before, with T = J r -M*. Mind that T is purely imaginary this time. 

Finally, in the unipotent case (type I with a = 0, or type III), normalize £ by setting for instance 
£(0) = i, £(0) = 0, so that £ is purely imaginary. The same function- and time-change yields \^\ - C, 

hence a natural basis of solutions is given by x\{9) = ^(9), X2{9) = ^(0) j Q To get globally define 
solutions, one avoids the double zeros of £ by drawing half-circles in the upper half-plane. Then the 



monodromy matrix is f J ^ J, with T = Jf, ^ = Jf, 



3.3 Spectral decomposition of the model operators 

We shall need below the spectral decomposition of the three model operators -\(d 2 x -a 2 x 2 ), -\{d 2 +a 2 x 2 ), 
-\d 2 x introduced above. They are essentially self-adjoint on C^°(R) by the classical Sears theorem (see 
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131 , Theorem 1.1 chap. 2 for instance), so the spectral theorem applies. The first operator has a pure 
point spectrum, while the second and the third have an absolutely continuous spectrum. Note that - jd 2 
is non-negative, while the spectrum of -\(d 2 + a 2 x 2 ) is the whole real line, as the following Lemma 
proves. 

Lemma 3.3.1 



1. (elliptic case) 



The spectral decomposition of L 2 (W)for the operator -\{d 2 x - a 2 x 2 ) is given by 



L 2 (R) = W-^ij (3-25) 

where L 2 ^ + , ^ is one-dimensional, generated by the normalized Hermite functions Ca l ^ 4 e~ ax2 ^ 2 He n (x 
for some constant C (see for the notations and normalization). 

2. (hyperbolic case) 
Set, for A e R, 

1/4 



12 

mix) : 



e XI%a e -\a^H 

1 1 iA 1 7 2 V« 1 i/1 3 7 

i^i (t(1 + — ), ~nax 2 ) ± , .. e 17r/4 x . iFi(-(3 + — ), -;iax 2 ) 



[r(| + ^) i a- 2' ' + - 1 "4 V a" 2 

where \F\ is the usual confluent hypergeometric function. Then Hft* = At]/*, and the (ft*, A e R) 
form a complete orthonormal system of generalized eigenfunctions of the operator H - —\(d 2 + 
a x ), so that any function f e L 2 (R) decomposes uniquely as 



f(x)= I ft + A (x)g + (A)dA + I m(x)g-(A)dA (3.27) 



with g ± (A) = ^f(x)ft*(x) dx. In particular, the following Parseval identity holds, 

f \f(x)\ 2 dx= f \g\A)\ 2 dA + f \g-(A)\ 2 dA. (3.28) 

Jr Jr Jr 

3. ( unipotent case ) 

Set, for A > 0, ft*(x) = e ±lxy ^_ Then Hi]/* = Aft* and the ft*, A > 0,form a complete orthonormal 
system of generalized eigenfunctions of the operator H = -Id 2 , with the usual Parseval-Bessel 



identity. 



Proof. 



(i) is classical and (iii) is straightforward by Fourier inversion and the usual Parseval-Bessel identity. 
Case (ii) is less common, though it can certainly be found somewhere in the literature. Let us explain 
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briefly how to obtain its spectral decomposition for a = 1. The easiest way is to remark that H = AAA 

„ „ . ., „ i . , , ■ , • • • / cos7r/4 -sin/r/4 

where A = k{xo x + o x x) = i(xo x + i) and A is the image ol the rotation matrix 1 



-l 



sin7r/4 cos7r/4 / 

by the metaplectic representation. The operator A is unitary. Explicit formulas found for instance in lfT31 
show that 

(Af)(x) -iV2e i7r/4 e i7rv2 e~' m{x ^~ y) 2 f(y) dy. (3.29) 
Jo 

As for the operator A, it is conjugate to i(d y + ^) after the obvious change of variable x = ±e y , hence 
its spectral decomposition is given by Fourier inversion on either half-lines, A(p* = A(p* (A € R) with 

(p*(x) = x ± 2 constituting an orthonormal basis of generalized eigenf unctions. Finally, t/f* := Acf>* may 
be obtained by applying the following formula (see lPT0ll ) 

x v - l e~^- yx dx = (2 y 6)- v/2 T(v)^ 2/8 ^D„ v (^) 

Jo 

(Re /3, Re v > 0) where D v is a parabolic cylinder function, also given by 

DAz) = T»e** i-T^- i;z 2 /2) - z-^- ,F^, ^/2)1 (3.30) 

[T(i(l-v)) 2 2 T(-v/2) 2 2 J 

(see [7], 8.2. (4) p. 117). □ 

3.4 Monodromy of non-resonant harmonic oscillators (elliptic case) 

We assume here that D e S°ff * s °f class (i) with a > 0. Then D is conjugate by a transformation in SV 
to an operator of the type -2ide - d 2 + a 2 x 2 + y where a > and y is a constant. Choose £ = - so that 
satisfies a normalized Pinney-Milne equation. Then Proposition 3.1.4 shows the following: 

Theorem 3.4.1 

The solution of the Schrddinger equation with arbitrary initial state 

<A(0) := J^CnhniO) (3.31) 

n>0 

is given by 

t A(0):=2]c^- i (' !+ 5)^/2/ J , 1 (0). (3.32) 

n>0 

The monodromy operator is given by the 'infinite-dimensional' monodromy matrix Md :— diag(e lA ",n € 
N), with A„ := -2n(n + j)a - ny. 

3.5 Monodromy of harmonic repulsors (hyperbolic type) 



One assumes now that D e S^J is either of class (i) with a < or of class (ii). Consider again the 
Ermakov-Lewis invariant 



x 2 + (tfd x + ^x) 2 



(3.33) 
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where one has assumed that £ = \r\ is purely imaginary this time, and /y 2 (£) - 2. Note that &L{iif) is 
anti-hermitian. Then 



££(177) ~ i* _ _1 
i?7 2 



„ 2 • *7 o ^ 1 ~ W 2 1 • V ^ 2fc 

- l-xd* + 5 — x - -l- + — 

?7 77 2 77 77 



(3.34) 



Suppose ^ ^ is an eigenvector of the Ermakov-Lewis operator, £>£,(irf)if/k - ikif/k- Then Propo- 

15- 

"4 7/- 



sition 2.1.4 implies that 1/^ := exp-4^x 2 . 1/^ is a generalized eigenfunction of the model harmonic 



repulsor, namely 



1 x 2 k 

2 77 z n 



(3.35) 



Hence: 
Lemma 3.5.1 



1. The equation (S£,{in) - ik)if/k - (k e R) /ms ftvo linearly independent solutions, 

tf ven (e, x) = y/TQbif'We&e-** . ^(1(1 + 2i*), I; ^) 

r(i + i|) 4 2 77 

- 2 V^)i/V/4 e il*V^ 1^ - ^(1(3 + 21*), |; — ). 

r(| + i|) 4 2 77 

The functions {{ilj k even ,ijj k odd ),k e R) constitute a complete orthonormal system for the operator 
£X(i77). 

2. f?7ie /las 

fl77<i 

Hence x — > -^exp|/t J" 6 * j^venW aw< ^ x — * |exp^/c J° ^rjj^ k odd (x) are solutions of the 
Schrodinger equation. 

Proof. 

1. is a direct application of Lemma 3.3.1, while 2. follows from an easy computation using the 
confluent hypergeometric differential equation zj^j iFi(a,c;z) + (c-z) d \ \F\(a,c,;z)-a \F\(a,c;z) = 0. 
□ 

The eigenf unctions i// k ven , ty k odd depend analytically on ^ for <f € C \ R_. If the operator D is of type I 
(so that £ has no zero), say with 7 = 0, then the phase exp ( k J -j?^ J gives the monodromy. If D is of type 
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II, then one must resort to a deformation of contour in order to avoid the singularities, as in the classical 
case, see section 3.2. Mind that the deformation of contour may change drastically the behaviour of the 
functions i// k ven , ty k odd for large x or large k (for instance, i^ ven and ip k dd become exponentially increasing 
for large x). Hence, in order to be able to follow the phase shift of the eigenfunctions ^ k even , if/ k odd along 
the contour T without getting divergent integrals, it is better to assume to begin with that the 'Fourier 
transform' (with respect to the spectral decomposition of £X(<f)) of the solution has compact support. In 
other words, the solution of the Schrodinger equation with initial state 

iA(0,x):= \ c + (k)if/ c even (0,x)dk + \ c-(k)i/f odd (0,x)dk 
Jr Jr 

for z e r (complex time), where c+, c_ are assumed to be compactly supported, is given by 

X) = C + (k)e k % ^Yeveniz, *) dk + || jT C-(k)e^ ^^(z, x) dk 

An immediate corollary is: 
Theorem 3.5.2 

Let i//(0) € L 2 (R), with decomposition 

<A(0,x):= f c + (k)tf ven (0,x)dk + f c-(k)i/f odd (0,x)dk. (3.36) 
Jr Jr 

Then the solution of any type (ii) Schrodinger equation with initial state if/(0) is given at time 9 = 2n by 

W2x,x)= f c + (k)e kT -^ even (0,x)dk+ f c-(k)e kT -^ odd (0,x)dk (3.37) 
Jr Jr 

where T = ^ or J r ^ (depending on the class of Vj), with T chosen as in section 3.2, is purely 
imaginary. The associated monodromy operator in S(L 2 (R), L 2 (R)) is unitarily equivalent to the multi- 
plication by the function k — > e kT ~ m7 with modulus one. 



3.6 Monodromy of non-resonant operators of unipotent type 

Suppose now D e S a ^ is of class (i), a = or (iii). Then 



(%d x + y x f 



(3.38) 



(£ e Staby 2 ) is an invariant of D (note the difference with respect to Proposition 3.1.4). Case (i), a = is 
trivial, for it is conjugate to the free Schrodinger equation. So assume D - -2idg - d 2 x + V2X 2 is of class 
(iii). Take £ = ir/ with 77 > as in section 3.2. Then (if k > 0) 



= —z{d x ---x) --. 

f; 2 2 77 77 



(3.39) 
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So 

i f] 2 1 2k 

if'k±(x) '•- exp ——x . exp±i, — x (3.40) 
4?7 ^77 

constitute a complete orthonormal system for £X(i?/) (the same statement holds true for potentials of 
class (i), in which case rj = 1 and the exponential prefactor is trivial). A short computation shows that 



2k i rj\ 
Di//k,± = ^~WK±- 



Hence one has the following: 
Theorem 3.6.1 

Let ip(G) € L 2 (R), with decomposition 



iA(0,x):= f c + {k)if/ h+ {x) dk + [ c-(k)if/ k -(x) dk. (3.41) 

Jr+ Jr+ 

Then the solution of any type (i), a = or type (Hi) Schrddinger equation with initial state i//(0) is given 
at time 9 = 2n by 



n,x)= f c+(k)e kT - iny ip h+ (x)dk+ f c-(k)e kT -' wy il/ h -(x)dk, (3.42) 



where T = ^ or J Y Y chosen as in section 3.2 (depending on the class of V2) is purely imag- 
inary. The associated monodromy operator in S(L 2 (R), L 2 (R)) is unitarily equivalent to the unitary 
operator on L 2 (R + ) given by the multiplication by the function k — > e kT ~ my . 



4 Symplectic structures and general solution of the Schrodinger equation 

The general emphasis in this section is so to speak on the non-quadratic part of the potential, namely, on 
Vq and V\ if D = -2ide - d 2 x + V2(9)x 2 + V\(9)x + Vq(9). It contains somewhat loosely related results: a 
definition of a three-dimensional invariant (£, 5\, 62); a generalization of the Ermakov-Lewis invariants to 
general potentials; a symplectic structure on a space 'containing' sucn that the SV-action becomes 
naturally Hamiltonian; finally, the computation of the monodromy for the 'resonant' operators of type 
(i)bis, (iii)bis. 

Definition 4.1 

We shall say that D e S a / 2 ^ is of generic type if: D is of class (i), D conjugate to D a ^ y = -2idg -d\ + 
ax 2 +y with a + n 2 /4, n - 0, 1, . . .; orD is of class (ii), D conjugate to D„ iQ , y - -2ide-d 2 + u n ^ a (9)x 2 +y. 

Denote by <S^f/ gen the set of operators of generic type; it is a disjoint union ofS V -orbits. 

Note (see Corollary 2.4.4) that the isotropy group of an operator D of generic type is generated by 
M\ and some £^ + }) ' f y + Mf 2 with £ + 0. 

Definition 4.2 
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Let D = -7id e - d\ 
Define: 



+ V 2 (9)x 2 + Vi(0)x + Vb(fl) 6 S[ 



fff 

<2,gen 



be of generic type. 



( i) £(D) to be the unique (up to a sign) periodic vector field such that £(D) e Staby 2 and 7y 2 (£(D)) = 2 
(£ real in the elliptic case, purely imaginary in the hyperbolic case); 

(ii) 8\(D) to be the unique periodic function such that 



Observe that -£^ + J/^ + Mg 2 e Lie(Go) is indeed unique (up to the addition of a constant times Ati) 
as follows from Corollary 2.4.4. The ambiguity in the definition of 82 may be solved by choosing for each 
SV-orbit an arbitrary base-point, an invariant (£81,82) for this base-point, and transforming (£,61,62) 
covariantly by the adjoint action along the orbit. Some non-local formulas fixing 62 more explicitly can 
probably be found, at least for potentials of type (i) (see Lemma 2.4.5), but we shall not need them. 

Another problem comes from the fact that the map (V2, V\, Vq) — > (£81,82) is not one-to-one (nor 
onto). Suppose one has some triple of functions (£81,82). Under some conditions that we shall not 
write explicitly (depending on the class of the potential), (£,61,62) is an invariant for some potential 
(V2, Vi, Vq); the quadratic part V% is given (by definition) by V2 = 2^(2 - ££ + \i; 2 ). (Supposing £ has 
only a finite number of zeros, all of which are simple or double, one has some rather straightforward 
conditions on the values of g| and -j* at the zeros of £ that ensure that £ e Staby 2 for some potential V^)- 
But V\ is not determined uniquely if £ does not vanish on the torus, since £ -3 ^ 2 is in the kernel of the 
operator £d + l£ (see formula (14.11) ). This can easily be explained by supposing (by conjugating by some 
element g e S V) that D is the model operator D = -2idg-d x +ax +y (a generic). Then £ is proportional 
to the constant vector field Xi which commutes with J/i, hence the invariant (£,61,62) is left unchanged 
by space-translations, whereas the operator D (and also the generalized Ermakov-Lewis invariant defined 
in Theorem 4.4 below) is not. Hence the vector invariant (g, 61,62) parametrizes Schrodinger operators 
of type (i) 'up to space-translations'. On the other hand, the map (V2, V\, Vq) — > (£ 61,62) is one-to-one 
for operators of type (ii) (up to a sign for £). 

It is not a priori self-evident that 82 defined by equation (14.21 ) is a periodic function. Considering the 
'inverse problem', i.e. supposing that the invariant (£61,62) is given, and supposing £ does not vanish 
on the torus, one must also check that every choice for V\ gives a function 62 which is periodic. This is 
the content of the following lemma: 



6\(D) + V 2 6i(D) = ~(Vi£(D) + \vit(D)); 



(4.1) 



C Hi) 



62(D) to be the unique periodic function (up to a constant) such that 




(4.2) 



Lemma 4.3 



One has: 



(4.3) 
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This formula implies: J^" £ 3 ^ 2 6i - 0; J^ n V\6i - 0. 
Proof. 

Using the invariant equations & - ^ 2 + 2V2$ 2 - 2 and 8\ + V261 = -\(V\% + \V\^), one obtains 

d 1 . 1 d 1 . 1 d ~>n in 



hence the first equation, which implies immediately: J^ n g 3/2 (8)di(6) d9 - 0. Hence (considering the 

inverse problem), if some potential V\ verifies J Vi(#)(5i(#) d6 = (so that 52 is well-defined), then 
this is also true for all possible potentials V\ . Now, integrating the first equation, one gets 



hence 



Hence 



£ 5 Vi - -2 



vtfwimdtf = J (^Vi)(e')(r 5 «5i)(0')d0' 

- -[(rwi(0)) 2 + (J r 3/2 (^)<5i(^) ^ 

and the integral over a period is zero. 

The following covariance result is an extension of Theorem 3.1.6. 
Theorem 4.4 



(4.4) 



Let D e S^/ gen be of generic type, with associated invariant (£ = %(D),5\ - 5 {(D), 62 - 62(D)). 
Then: 



&£(D) := - 



1(1 + I^ )x 2 _ ^2 + ^ {xdx + QxX) + ^26 { (-id x ) + (Vtf + 26 1 )x) + 2(5 2 + \v<£) 



(4.5) 



is an invariant for the Schrddinger operator D. 

2. Let ((p; (a, b)) e SV and g : (6, x) — > (9', x') = ((f)(6), x ^(p(6) - a(6)) be the associated coordinate 
change. Then 

n llA (<P; (a, b))&£(D)7T U4 (4>; (a, b))~ l = &£(D) (4.6) 

where S£(D) is obtained by applying the transformation g to the coordinates, changing the po- 
tentials Vq and V\ by the cri^-action ofSV, and transforming the invariant as follows: 



(4.7) 
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Sx = 0* . 8y o tfr 1 + (fa - -al)\ (4.8) 

62 = 62 o (fT l + {8\d - aSi) + ffc + (|(a 2 - ad) - \aa - ga 2 ). (4.9) 
Furthermore, (g, 8\ , 82) is the invariant associated to tr\ /4(D). 

Proof. 

1 . Look for an invariant of the form 

- [a(9)x 2 - b(9)d 2 x - ic(G)(xd x + d x x) + d(6)(-id x ) + e(6)x + f(9)] (4. 10) 

and solve in a, b, c, d, e, f. One obtains the following constraints: 

d = 2V 2 c, b = -2c, c = -a + V 2 b (4.11) 

- whose general solution is in Proposition 3.1.4 above, namely, a = |(1 + \^ 2 ), b - c - —jg - 
and the set of following equations: 

d = Vib-e, e = V l c + V 2 d, f = -dV\ (4.12) 
which implies the compatibility condition 

2. Since (assuming = 2 is fixed) there is a unique invariant for operators of generic type, one 
necessarily has 

£| + V gi + Mg 2 = Ad{<P; {a, b)).{£ ? + & Sl + M Sz ) (4. 13) 

which gives the above formulas for (f, 81,82). 

There remains to check for equation ( 14.6I ). Consider first the covariance under a time-reparametrization 
(p. It has already been proved for the quadratic part of the Ermakov-Lewis operator, see Theorem 
3.1.6. The linear part -2(-i8id x + (V\% - 8\)x) transforms covariantly under <f> since (see proof of 
Theorem 3.1.6) 

V ^~^= <P'~kVi.^-8 { - i£si), -\S x d* = -i5 1 d x (4.14) 
d9 2<p 

and 

= ^^ e -4 x \-i 6l d x + -tfOx^V'^Wfl.*^^ 

= (-iStf* + (Vi| - SirfWP, A (4-15) 
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As for the zero-order term -M62 + \Vq£), it is obviously invariant under the conjugate action of 
7Ti/4(0). Since Vol - (Vof) <A _1 > this implies also 62 = 82 ° <f>~ 1 . 

Consider now the covariance under an infinitesimal nilpotent transformation J/* + M/ 2 . One has 
[ad x + iax + b, ££(£, S\ , 6 2 )] = 

j[ad x + iax, 1(1 + -^V - + i^xd x - 2(6 1 (-id x ) + (Vrf - 

j (1 + + of) x - {at + 2at)(-id x ) - 2a(Vi{ - <5i) + 2as\ , (4.16) 

to be compared with the infinitesimal change of &£, under the transformation x — > x + sa, 61 — > 
5i + e(^a - 5«|), £2 — > #2 + e(((5ifl - a#i) + %b), V\ — > V\ - 2s(a + V2a). This is a straightforward 
computation, which requires the use of the equation defining namely, £ = 4(1 + |^ 2 ) - 2V2%- 

□ 




Using the parametrization of en by the vector invariant (£, <5i, #2), one can easily define a natural 
symplectic structure on a linear space Q. and a hamiltonian action of S V on f2 reproducing the SV-action 
on ojr? . 

<2,gen 

Definition 4.5 

Let Q, - C°°(K./27rZ, R 4 ) be the linear manifold consisting of all 2n-periodic vector-valued C°° func- 
tions X(t) := (p, q, E, t)(r), r e R/27rZ with singular Poisson structure defined by 

{ P (t), q(T')\ = 5{r - t'), {E(t), t(T')\ = Sir - r') (4.17) 



See for instance lfl4l . chap. X for some remarks on distribution-valued singular Poisson structures 
on infinite-dimensional spaces. The energy E is canonically conjugate to t, which allows us to consider 
generalized canonical transformations for which t is a coordinate. This usual trick for hamiltonian sys- 
tems with time-dependent Hamiltonians can for instance be found in |9). Hamitonian vector fields Xe, 
for H = H(p, q, E, t), acts separately on each fiber t =constant, namely, 

(X H f)(T) := {{dpHdg - d q Hd p + d E Hd p - d t Hd E )f] (r). (4.18) 

Definition 4.6 

Let (£, 5 1,82) be a triple of 2n-periodic functions. Define O := (D(£, <5i, 82) to be the following 
functional on fl, 

<®, X) = <£ L(t(f))E(r) + ^(t(T))p(r)q(r) + Si (Kt)) P (t) - 6, {t(r))q{r) + S 2 (t(r))\ dr. (4.19) 
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Theorem 4.7 

Represent £,f + 3f g + Mh £ so by the hamiltonian vector field Xfi{f,g,h) associated to 

H(f,g, h) := -(f(t)E + X -f(t)pq + -Jq 1 ) - (g(t)p + g(t)q) - h(t). (4.20) 

Then the action of Xh on the functional <!>(£, 5\, 62) coincides with that given in Theorem 4.4. 
Proof. 

Observe that the map from st> to the Lie algebra of vector fields on Q given by X/ + }f g + Mh — > 
XH(f,g,h) is a Lie algebra homomorphism. The vector field Xh is given explicitly by 



XH(f,g,h) 



l -f{t){qd q - pd p ) + f(t)d t - X ~f(t)qd p 



+ 



i\f\t)pq + fm + X -f"\t)q 2 ) + (gp + gq) + h(t) 



[g(t)d q ~ g(t)d P ] 
d E - 



(4.21) 



The rest is a straightforward computation. 



Let us conclude this section by computing the monodromy for 'resonant' operators of type (i)bis and 
(iii)bis. 

Consider any resonant operator D. The associated classical monodromy is unipotent. We choose 
£ e Staby, to be purely imaginary, £ := irj as before (see section 3.2). A generalized Ermarkov-Lewis 
invariant may then be defined as 



(i,*d, + Uxf 



+ - [d(-id x ) + ex + f], 



(4.22) 



where d, e, f are defined as in Theorem 4.4 but with £ replaced by 77 (see equation (14.101) for notations). 
Hence 



6£(D) - ik _ 1 



i fi 9 d e f 

2r/ t] T] t] 



k 

n' 



(4.23) 



Suppose &£,(D)if/k - ikftk an d set 



, i V o id 
iffk = exp\---x +--x\if/ k . 
4n 2 77 



(4.24) 
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Then a simple calculation gives 



1 , f] 

-d— + - \x + 

2 rj z rj 



-f + 2k lid 



(4.25) 



If D is of type (i)bis, then d, e, f (easy to obtain from Theorem 4.4 and the isotropy algebra given 
in section 2) satisfy \d\ + | = identically, so the model operator is (up to a constant) the Laplacian 
as for case (iii). Then the monodromy can be computed along the same lines as in section 3.6, with a 
time-independent shift in k due to the function -/ + ~ - "iff^ 1- 

Lemma 4.8 

Let D = -2idg - d 2 + n 2 x 2 + C cos n{9 - cr/2)x + y be a Schrddinger operator of type ( i)bis. Set 



with d = -£ sin 3n(6 - cr/2), 77 = 1 - cos(2«6» - cr), k' = k + 3 (-fi-f- Then 

. . \ 

j -3/2 ,/^77 1 *7 



£><A*,± ^ 



'2*: 














27, 



■y 



(4.26) 



(4.27) 



Proof. Tedious computations. □ 

Apart from the time-periodic shift j (fy = (-j^j) X^e-oV^ ( wn i cn * s integrable on the contour T) 
and the time-independent shift in k, one is left once again with a phase proportional to k/rj (note that the 
term in drj~ 3/2 ^/2~F is irrelevant since J Q (dr]~ 3/2 )(6) d9 = by Lemma 4.3; recall d - -26 1 by Theorem 
4.4). 

Hence one obtains: 
Theorem 4.9 

Let i//(0) € L 2 (R), vwY/z decomposition 

i//(0,x):= J c + {k)i// K+ {0,x)dk + J c-(k)i// k -(^ x) dk. (4.28) 

Jr+ Jr+ 

r/je« solution of the type (i)bis Schrddinger equation 

{-2id e + d 2 x + n 2 x 2 + C cos(n9 - cr/2).x + y)\j/ = (4.29) 
with initial state tA(0) w given at time 6 = 2n by 

ijj(2n,x)= f c + (k)e k ' T - i ' Tf ^ K+ (0,x)dk+ f c-(k)e k " r - brf ^ K -(0,x)dk, (4.30) 

Jr+ Jr+ 

where k' - k + 3 ( jfjj) » ^ = J) 2 " is purely imaginary) and 

7 = 7+7 f(~) ( e ) de - 
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The associated monodromy operator in S(L 2 (W), L 2 (K.)) is unitarily equivalent to the unitary operator 
on L 2 (K.+) given by the multiplication by the function k — > e kT ~ m ^. 



Suppose now D is of type (iii)bis. Then the x-coefficient in the transformed Ermakov-Lewis operator 
(14.251) does not vanish, so one must take for 'model operator' -d 2 + x, whose eigenfunctions are related 
to the Airy function. The solution of the monodromy problem will be given by a series of lemmas. 
In the sequel, 77 = (1 + sin«#)(l + asinnO) is the (real- valued and non-negative) invariant, and r] 1 ^ 2 = 
(1 + a sin nG) 1 ^ 2 cos(| - n|)x is the smooth square-root of 77 chosen in section 2. 

Lemma 4.10 

Let Ai be the entire function, solution of the Airy differential equation (—d 2 + x)Ai(x) = 0, defined on 
the real line as 

Ai(x) 



- f cosf— + xt\ dt. (4.31) 
Jo \ 3 / 



It is (up to a constant) the only solution of the Airy differential equation which do not increase exponen- 
tially on R+. The functions fk(x) := Ai(x — k), k el define (up to a coefficient) a complete orthonormal 
system of generalized eigenfunctions of the self-adjoint closure of the Airy operator —d 2 + x with core 
C™(R) c L 2 (R). 

Proof. Easy by using a Fourier transform. □ 
Lemma 4.11 

The x-coefficient in the transformed Ermakov-Lewis invariant (14. 251 ) reads 

LI + f = -C aV ^ 2 (4.32) 
2 r\ l n 



where C a = (1 - or)(l + a 12) Vl - a 2 . 
Proof. 

Computations similar to that of Lemma 4.3 (with the simple difference that $ - ^ 2 + 2V2^ 2 = 
here) yield 

d 1 C n 

-(77-^) = 77 - -2 (4.33) 
ad 77 

where C a is some constant which must be chosen in order that the right-hand side be 27r-periodic. Note 
that the singularities in the above equation are only apparent; one may avoid them altogether by using a 
contour T in the upper-half plane as in section 3.2. Since f 77 = 27r(l + a/2) and f -7^ = 2n , — = 

(see Proposition 2.2.3), this gives C a — (1 - a)(\ + a/2) - a 2 . Then a straightforward computation 
yields formula (14.321) . □ 

Lemma 4.12 

Set _ 2 

MO, *) = exp [^-x 2 - ~xj .17-! A/(xCyy ^ - C~ 2l \-f + 2k + I— )J (4.34) 
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Then 



2k 



+ 



Li lit 

2t] 2 U 



(4.35) 



Proof. 

The <A/t are obtained as in section 3.5 (monodromy of hyperbolic operators) by taking a complete or- 
thonormal system of generalized eigenfunctions ^ for the transformed Ermakov-Lewis invariant (14.251) 
and going back to the functions i/^. Then (14.351) is proved by a direct tedious computation. □ 

One may now conclude: 

Lemma 4.13 

Let i/f(0) € L 2 (R), with decomposition 



i/f(0,x):= J c(k)if/ k (0,x)dk. 
Then the solution of the type (iii)bis Schrddinger equation 



(l 7t \ 

-2ide + d\ + v„ jQ .x 2 + C(l + asinnO)? cos(- - n-).x + yj\p - 



with initial state tA(0) is given at time 6 = 2n by 

ijj(2n, x) = J c(k)e kT - in ^M0, x) dk, 
where T = 4^- (T is purely imaginary) and 



y = y + 



I(-H(; 



(9) dO. 



(4.36) 



(4.37) 



(4.38) 



The associated monodromy operator in S(L 2 (K), L 2 (R)) is unitarily equivalent to the unitary operator 
on L 2 (R) given by the multiplication by the function k — > e kT ~ m ^. 
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